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gain,  and  directivity  can  be  calculated. 

Computer  programs  based  on  this  new  formulation  were  used  to  deter¬ 
mine  the  impedance  versus  frequency  for  various  loop  antennas.  Verifi¬ 
cation  of  the  calculations  with  experimental  results  showed  that  the 
technique  is  highly  accurate  not  only  for  the  impedance  but  also  for 
the  current  distribution.  Calculated  and  measured  results  for  several 
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ANALYSIS  AND  DESIGN  OF  SPECIAL  ANTENNA  CONFIGURATIONS 


I.  OBJECTIVES 

The  objective  of  this  program  was  to  develop  digital  computer 
techniques  for  analyzing  the  Impedance,  efficiency  and  pattern  character¬ 
istics  of  wire  antennas.  In  particular,  studies  were  to  be  made  of  halo 
and  multi -turn  loop  antennas. 


II.  INTRODUCTION 

In  the  past,  antenna  development  has  progressed  via  experimentd i 
methods  guided  by  intuition  <.nd  experience.  An  antenna  is  built  and 
Its  characteristics  are  measured  in  the  laboratory.  The  antenna  dimen¬ 
sions  and/or  shape  are  then  modified  and  the  measurements  are  repeated 
until  satisfactory  performance  is  obtained.  This  approach  is  slow  and 
tedious.  Accurate  measurements  require  skill,  care  and  sophisticated 
instrumentation.  The  final  antenna  design  is  seldom  optimized,  since 
It  represents  simply  the  best  of  all  tha  models  tested. 

Now  the  experimental  approach  cannot  and  should  not  be  eliminated, 
but  It  should  be  auom .nted  and  coordinated  with  modern  computational 
techniques.  Indeed,  the  state  of  the  art  now  permits  computer  simula¬ 
tion  of  many  laboratory  experiments,  and  results  can  be  predicted  with 
greater  speed  and  accuracy  and  lower  costs.  With  a  suitable  program, 
a  modern  computer  can  deliver  a  design  in  one  day  that  might  require 
several  mor.chs  with  the  experimental  cut-and-try  approach. 

The  ElectroScience  Laboratory  has  developed  significant  capabil¬ 
ity  and  experience  in  the  analysis  of  thin-wire  antenna  problems. 

Thiele'* has  developed  and  published  a  solution  for  the  current  dir. tri- 
butions  and  radiation  patterns  of  a  Yagi  array,  and  is  now  able  to 
calculate  the  impedance  as  well.2  Richmond3  carried  out  a  computer 
analysis  of  an  array  of  V-antennas  for  a  NASA  satellite  system. 

Agrawan  measured  and  calculated  the  impedance  of  a  planar  array  of 
three  TEM-line  antenn.-r  interconnected  with  transmission  lines,  and 
obtained  excellent  results.  Jennetti  and  Agrawal6  measured  and  calcu¬ 
lated  patterns  and  impedance  of  a  broadside  array  of  TEM-line  loops. 

To  optimize  the  array,  they  programmed  a  computer  to  analyze  10,000 
different  designs  per  minute  and  display  the  best  results.  A  purely 
experimental  approach  would  have  required  several  months  to  achieve  the 
excellent  bandwidth,  VSWR  and  directivity  of  their  final  design. 

Thiele,  Travieso-Diaz  and  Jones6  measured  and  calculated  the 
radiation  patterns  of  a  dipole  antenna  mounted  on  the  base  of  a  conduct¬ 
ing  cone.  In  their  solution,  a  computer  is  programmed  to  enforce 
boundary  conditions  at  many  points  on  the  cone.  Tsai  and  Otto7»8 
employed  the  iourier  t>. '.inform,  discrete  boundary  matching  and  counter 
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optimization  to  calculate  the  impedance  of  dipoles  mounted  on  spheres 
and  cylinders.  Their  results  agree  closely  with  measurements. 

The  ElectroScience  Laboratory  has  made  significant  progress  in 
the  past  year  in  the  investigation  of  wire  antennas  for  the  Ballistic 
Research  Laboratory.  This  progress  is  surveyed  in  the  next  sections, 
and  the  details  may  be  found  in  the  Technical  Reports  included  as 
Appendices. 


III.  TASKS  1  AND  3:  MULTI -TURN  LOOPS 

A  new  reaction  solution  was  formulated  for  three-dimensional 
wire  antennas  located  in  free  space  or  over  a  ground  plane.  In  this 
approach,  a  piecewise-sinusoidal  expansion  is  employed  for  the  unknown 
current  distribution  on  the  antenna.  To  determine  the  expansion  co¬ 
efficients,  we  force  the  antenna  current  function  to  have  the  correct 
reactions  with  a  set  of  interior  test  sources.  This  yields  a  set  of 
simultaneous  linear  equations.  These  equations  are  solved  on  a  digital 
computer  to  determine  the  current  distribution,  and  it  is  found  that 
the  solution  is  stationary.  Then  it  is  not  difficult  to  calculate  the 
impedance,  radiation  efficiency,  field  patterns,  gain  and  directivity. 

Computer  programs  based  on  this  new  formulation  were  found  to  be 
exceptionally  fast,  accurate  and  versatile.  These  programs  were  used 
to  determine  the  impedance  versus  frequency  for  various  loop  antennas. 
Technical  Reports  2708-1  and  2708-5  (Appendices  A  and  B)  show  the  results 
for  several  one-turn  polygon  loops.  The  accuracy  of  the  calculations 
was  verified  by  comparison  with  experimental  measurements  for  the  square 
loop  and  with  Storer's  solution  for  the  circular  loop.  It  was  found 
that  our  new  technique  is  highly  accurate  not  only  for  the  impedance 
but  also  for  the  current  distribution. 

Report  2708-6  (Appendix  C)  presents  experimental  and  theoretical 
results  for  two- turn  and  three-turn  loops  over  a  ground  plane. 


IV.  TASKS  2  AND  4:  HALO  ANTENNAS 

In  developing  the  techniques  and  computer  programs  for  halo 
antennas,  preliminary  studies  were  necessary  with  simplified  configura¬ 
tions.  Technical  Reports  2708-2  and  2708-3  (Appendices  D  and  E)  show 
the  results  for  coplanar-skew  and  nonplanar-skew  cylindrical  wire 
arrays. 

Report  2708-4  (Appendix  F)  shows  calculated  and  measured  results 
for  several  three-dimensional  dipoles  and  a  square  halo  antenna.  The 
calculations  show  excellent  agreement  with  the  measurements.  This  report 
also  develops  in  detail  our  theoretical  formulation  for  three-dimensional 
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wire  antennas  with  finite  conductivity,  multiple  excitation  and  lumped 
loading.  The  report  includes  the  computer  program  and  the  theory  for 
the  current  distribution,  impedance,  radiation  efficiency  and  patterns. 


V.  RECOMMENDATIONS 

As  indicated  above,  we  have  developed  excellent  computer  tech¬ 
niques  for  wire  antennas  located  in  free  space  or  over  a  ground  plane. 
These  techniques  are  also  applicable  to  slot  and  dipole  antennas  mounted 
on  a  rocket  or  projectile  with  axial  symmetry.  For  example.  Figure  1 
compares  our  measurements  and  calculations  for  a  dipole  mounted  on  a 
cylindrical  vehicle.  The  same  computer  program,  based  on  the  reaction 
formulation,  handles  a  variety  of  configurations  including  biconical 
antennas,  slotted  cones  and  slotted  cylinders. 

We  recommend  an  extension  of  this  program  to  apply  these  tech¬ 
niques  to  a  configuration  of  direct  interest  to  the  Ballistic  Research 
Laboratory:  a  loop  or  halo  antenna  exciting  a  slot  in  a  conical  rocket 
or  projectile.  Our  experience  in  this  area  is  illustrated  in  References 
6  through  8  and  in  Figure  1.  The  theory  and  computational  techniques 
would  be  applied  to  optimize  the  impedance,  patterns,  bandwidth,  effi¬ 
ciency  and  voltage-breakdown  characteristics  of  the  projectile-mounted 
antenna. 
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Theoretical  and  experimental  admittance  of  dipole  on  cylindrical  vehicle. 
The  theoretical  results  are  based  on  the  reaction  formulation. 
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ABSTRACT 


A  brief  summary  of  the  application  of  the  reaction  concept  to 
calculate  mutual  coupling  between  thin  linear  dipoles  is  given. 
Using  reaction,  a  basis  is  developed  for  a  multi  segment  piecewise- 
sinusoidal  approach  to  the  calculation  of  mutual  coupling.  This 
technique  is  applied  to  compute  the  input  admittance  of  a  planar 
polygon  loop  antenna,  and  results  are  presented  which  show  good 
agreement  between  these  calculations  and  similar  theoretical  values 
for  a  circular  loop  of  equal  area. 


16 


INPUT  ADMITTANCE  OF  PLANAR 
POLYGON  ANTENNAS 


INTRODUCTION 

In  recent  years  greater  attention  has  been  directed  to  the 
accurate  calculation  of  antenna  impedance.  In  the  area  of  coupled 
linear  antennas,  much  of  this  work  has  been  based  upon  the  induced- 
emf  technique  first  used  by  Carter.!*  Others,  includinq  Storer?  and 
Tai,3  have  attacked  the  problem  using  variational  techniques  with 
significant  success.  A  third  approach,  which  can  be  obtained  through 
variational  considerate  r'-’s  and  which  is  reducible  to  the  induced-emf 
technique  as  one  case,  is  the  reaction  concept  introduced  by  Rumsey/ 

The  first  two  techniques  mentioned  above  assume  that  an  antenna 
is  represented  by  two  rol linear  thin  cylindrical  seqments.  Differences 
between  these  first  two  techniques  center  around  the  choice  of  the 
current  distribution  function  to  be  used.  The  induced-emf  approach, 
which  assumes  a  piecewise-sinusoidal  current  on  each  of  the  two 
antenna  arms  fails  for  the  full -wavelength  dipole  case,  while  the 
variational  approach  of  Tai  yields  finite  terminal  impedance  for 
all  lengths. 

The  reaction  concept,  like  the  other  techniques,  has  in  the  past 
been  concerned  mainly  with  the  two-segment  representation  of  a 
linear  antenna. 

Harrington,  Chang  and  King  ,  and  more  recently  Richmond/  have 
successfully  reformulated  the  coupling  problem.  Their  aoproach  is 
to  represent  a  dipole  by  several  collinear  segments  on  each  of  which 
the  current  distribution  is  described  by  some  fairly  simple  function. 
Application  of  the  reaction  principle  results  in  a  system  of  linear 
equations  relating  the  unknown  current  coefficients  to  the  voltage 
at  each  junction  between  segments.  The  set  of  equations  is  then  solved 
with  the  aid  of  the  digital  computer. 

The  form  of  the  technique  developed  by  Richmond  is  general  enough 
to  be  applied  in  a  straightforward  way  to  a  much  larger  class  of 
antennas  than  the  coupled  thin  linear  dipoles  with  which  the  aoproach 
was  initially  concerned.  The  purpose  of  this  communication  is  to 
describe  the  technique  briefly  and  to  show  its  application  to  the  case 
of  a  planar  polygon  antenna. 


THEORETICAL  DISCUSSION 

In  any  reference  to  the  use  of  reciprocity  or  reaction  to  compute 
coupling  between  antennas,  an  important  consideration  is  the  form 


References  are  Listed,  on  page  25 . 
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of  the  current  distribution  that  is  assumed  to  exist  on  the  antennas. 

One  is  at  liberty  to  approximate  the  true  current  distribution  in  any 
way  he  desires,  but  the  amount  of  work  one  must  then  do  to  obtain 
accurate  results  depends  greatly  on  his  particular  choice  of  an 
approximation.  This  discussion  will  consider  a  piecewise-sinusoidal 
approximation  for  reasons  other  than  just  its  obviously  close  re¬ 
lation  to  the  true  current  distribution  on  an  electric  line  source. 

The  use  of  reaction  to  calculate  coupling  between  two  cylindrical 
dipoles,  where  each  is  composed  of  two  segments  with  an  assumed 
piecewise-sinusoidal  current,  is  identical  to  the  induced-emf  technique. 
This  method  relates  the  unknown  short-circuit  terminal  current  I2 
Induced  on  antenna  2  by  voltage  V-|  impressed  at  the  terminals  of 
antenna  1  (with  unit  terminal  current)  by 


(1)  Z,2I2  =  V, 


where 


Z12  =  “  I2  -2  *  -1  ds  * 

and  the  integration  is  over  the  surface  of  antenna  2.  Here  (E^Hi) 
denote  the  free-space  fields  generated  by  a  linear  test  source  located 
along  the  axis  of  antenna  1.  Eq.  (2)  is  obtained  by  enforcing  re¬ 
ciprocity  between  the  test  source  and  antenna  2.  In  most  practical 
situations  the  surface  current  density  0?  on  antenna  2  can  be  adequately 
represented  by  a  line  current  along  trie  surface.  This  reduces 
Eq.  (2)  to  a  line  integral, 


(3) 
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Schelkunoff  and  Friis®  give  the  cylindrical  polar  components  of 
the  free-space  field  of  an  electric  line  source  on  the  z-axis  with 
sinusoidal  current  distribution  as 


(4) 
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where 

(6) 


/pr.  , 


V  is  the  potential  due  to  end  charqes,  and  the  primes  denote  derivatives 
with  respect  tc  z.  Fig.  1  illustrates  the  geometry  related  to  Eqs.  (4) 


Fig.  1.  Geometric  relations  for  an  electric  line 
source  on  the  z-axis  with  an  observation 
point  at  (p,z). 


and  (5).  If  an  antenna  is  represented  by  a  collection  of  these  line 
sources  joined  end  to  end  in  some  fashion,  and  if  the  current  is 
forced  to  vanish  at  any  free  ends,  then  the  partial  derivatives  of  V 
in  Eqs.  (4)  and  (5)  contribute  nothing  to  the  field  expressions  for  the 
whole  structure  and  may  be  neglected  providing  the  end  charges  are  neoliqible. 
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The  collection  of  line  segments  reoresentinq  one  or  more  antennas 
forms  an  N- port  network ,  where  N  represents  the  number  of  junctions 
which  have  a  segment  on  each  side  (i.e.,  no  open  end  can  be  a  junction). 
An  equivalent  viewpoint  is  that  the  collection  consists  of  N  overlapping 
dipole  and/or  Vee  antennas  on  each  of  which  the  current  vanishes  at 
either  end. 

Let  an  ordinary  dipole  be  a  special  case  of  a  Vee  antenna.  Also, 
let  a  test  Vee  antenna  be  located  on  axis  with  terminals  at  junction 
i_  and  with  arms  the  same  lengths  as  those  of  the  ith  true  surface-current 
Vee  antenna.  Then  the  coupling  between  the  ith  junction  any  any  other 
junction,  say  the  j™.  junction,  can  be  calculated  by  enforcing  re¬ 
ciprocity  between  the  test  Vee  antenna  on  axis  at  junction  i_  and  the 
true  surface-current  Vee  antenna  whose  terminals  are  located  at  the 
ith.  junction.  This  is  done  through  Eq.  (3)  which  for  every  possible 
pair  of  junctions  i_  and  j.  generates  a  matrix  of  coefficients  [Z-jj]. 

These  are  the  coefficients  of  a  system  of  linear  equations  each  having 
the  form 


(6)  l  Zjj  Ij  =  Vi  t  1  =  1,  2>  N  . 
j 

By  assuming  a  known  voltage  source,  usually  a  unit  voltage  source, 
at  a  junction  of  interest  and  by  requiring  short  circuits  at  all  other 
junctions,  the  unknown  current  coefficients  at  all  the  junctions  can 
be  determined.  The  current  coefficient  at  the  junction  with  the  voltage 
source  is  the  input  admittance  of  the  system  if  the  voltage  source 
is  a  unit  source.  The  input  impedance  is  easily  obtained  as  the 
reciprocal  of  the  input  admittance. 

Richmond  realized  several  advantages  by  using  the  piecewise- 
sinusoidal  current  description  on  the  test  Vee  antenna  and  on  the  over¬ 
lapping  true  surface-current  antennas.  First,  the  reaction  between 
any  two  (nonintersecting)  coplanar  dipole  or  Vee  antennas  may  be  cal¬ 
culated  in  closed  form  in  terms  of  sine  and  cosine  inteqrals.  The 
arms  of  either  antenna  can  be  of  unequal  length,  and  their  Vee  angles 
need  not  be  equal.  Also,  a  good  value  for  the  mutual  couplinq  may  be 
obtained  for  individual  antenna  arms  up  to  a  quarter- wavelength  in 
extent.  Of  course,  greater  accuracy  may  be  had  by  using  more  segments, 
especially  when  the  test  source  and  the  true  surface-current  antenna 
are  relatively  close  together.  However,  this  technique  has  been  used, 
for  example,  to  calculate  the  input  impedance  of  a  linear  dipole  to 
the  same  degree  of  accuracy  but  with  far  fewer  segments  than  were 
needed  for  a  similar  technique  employing  a  different  current  distribution.7 

These,  then,are  the  advantages  of  the  piecewise-sinusoidal  multi¬ 
segment  approach  to  couplinq  between  coplanar  linear  antennas:  relative 
simplicity  of  concepts;  speed  and  accuracy  of  numerical  calculations; 
and,  as  the  next  section  attempts  to  demonstrate,  applicability  to  a 
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more  general  class  of  wire  antennas. 


PRESENTATION  OF  RESULTS 

The  technique  discussed  in  the  previous  section  has  been  employed 
successfully  to  calculate  the  self-impedance  of  linear  dipoles,  to  ob¬ 
tain  the  mutual  impedance  between  two  linear  dipoles  in  various 
orientations,  and  to  calculate  driving-point  impedance  and  current 
distributions  for  coupled  coplanar  Vee  antennas.7 

The  results  to  be  reported  here  deal  with  the  calculation  of 
driving -point  admittances  and  current  distributions  for  square  loops 
and  polygons.  Part  of  the  motivation  for  this  investigation  was  the 
desire  to  find  some  representation  for  a  circular  loop  antenna  that 
could  be  analyzed  by  the  reaction  technique  using  straight  wire 
segments.  The  regular  polygon  antenna  was  proposed  to  fill  this  need. 

The  criterion  chosen  to  relate  the  approximating  polygon  to  a  circular 
loop  was  that  of  equal  enclosed  areas.  There  are  at  least  two  reasons 
for  enforcing  this  relationship.  First,  a  circle  is  more  closely 
approximated  geometrically  by  a  concentric  N-sided  polygon  with  equal 
enclosed  area  than  by  either  an  inscribed  or  circumscribed  N-sided 
polygon.  Second,  the  dependence  of  both  the  field  expressions  and 
the  input,  impedance  of  the  circular  loop  antenna  on  its  enclosed  area 
suggested  that  this  feature  be  incorporated  in  any  approximating 
structure. 

Figure  2  is  a  plot  of  the  input  conductance  and  susceptance  of 
a  square  loop  antenna  as  a  function  of  frequency.  These  calculations 
bear  a  strong  resemblance  to  the  comparison  shown  in  Fiq.  3.  These 
plots  compare  the  input  conductance  and  susceptance  of  an  octagonal 
loop  antenna  to  those  of  a  circular  loop  antenna  with  equal  enclosed 
area.  The  circular  loop  data  were  taken  from  a  paper  by  Storer9  in 
which  the  input  admittance  was  obtained  tnrough  consideration  of  a 
Fourier  series  solution  obtained  by  Halle'n.  Notice  that  the  comparison 
extends  beyond  a  fuV.  wavelength  in  circumference. 

The  input  admittance  data  for  the  octagon  were  calculated  by  using 
two  segments  to  represent  each  side  of  the  figure.  There  was  negligible 
difference  between  admittances  calculated  by  placing  the  input  terminals 
at  a  vertex  of  the  octagon  and  those  calculated  with  the  input  at  the 
middle  of  one  of  the  sides.  In  addition,  the  agreement  shown  in  Fig.  3 
over  the  entire  range  of  frequencies  is  only  slightly  different  from 
that  obtained  by  using  an  approximating  octagon  havinq  only  one  segment 
per  side  with  input  terminals  at  a  vertex.  The  data  for  the  square 
loop  were  calculated  using  two  segments  per  side  at  the  higher  frequencies 
and  only  one  segment  per  side  at  the  lower  frequencies. 


.21 


CONDUCTANCE  (  mmhos  )  SUSCEPTANCE  (mmho*  ) 


Fig.  2.  Input  conductance  and  susceptance  of  a 
square  loop  vs.  circumference  of  a 
circle  (in  wavelengths )havinq  the  same 
enclosed  area.  R  is  the  radius  of  the 
circle,  and  k  =  2n/A. 
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-  CIRCULAR  LOOP  (STORER) 

•  •  •  •  16  SEGMENT  OCTAGON 


The  current  distribution  on  the  octagon  as  determined  from  the 
values  of  the  current  coefficients  at  the  junctions  between  segments 
shows  good  agreement  with  circular  loop  data  in  both  form  and  phase, 
but  the  magnitude  does  not  show  good  agreement  for  circumferences 
under  half  a  wavelength.  At  approximately  one  wavelength  in  cir¬ 
cumference,  however,  the  current  coefficients  at  the  junctions  of  the 
sixteen-segment  octagon  show  excellent  agreement  in  both  magnitude  and 
phase  with  corresponding  values  of  the  current  distribution  on  a 
circular  loop  antenna  as  published  by  Storer.9 


CONCLUSIONS 

The  multi  segment  reaction  technique  developed  by  Richmond  for 
coupled  straight-wire  antennas  and  coupled  Vee  dipoles  has  been  utilized 
to  calculate  the  input  admittance  and  current  distribution  of  the 
polygon  loop  antenna.  Moreover,  the  polygon  loop  has  been  shown  to 
approximate  a  circular  loop  for  the  puroose  of  calculating  input 
admittance  and  current  distribution. 

This  method  forms  the  basis  for  computing  the  impedance,  current 
distribution,  and  far-field  patterns  of  an  antenna  or  a  system  of 
antennas  composed  of  loops  or  of  combinations  of  loops,  such  as 
multiturn  loops,  and,  in  general,  any  antenna  made  up  of  or  approximated 
by  piecewise  straight-wire  segments. 
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ABSTRACT 


The  reaction  concept  is  applied  to  wire  antennas  in  the  shape  of 
regular  polygons  to  calculate  the  input  impedance  (or  admittance)  and 
current  distribution.  The  calculated  results  are  found  to  give  excellent 
agreement  with  experimental  measurements  and  independent  theoretical 
data. 


28 


COMPUTER  ANALYSIS  OF  POLYGON  LOOP  ANTENNAS 


I.  INTRODUCTION 

This  report  deals  with  a  class  of  wire  antennas  that  has  hereto¬ 
fore  been  largely  ignored.  Except  for  some  preliminary  results1  this 
report  constitutes  the  first  presentation  of  a  rapid,  efficient  technique 
for  calculating  current  distribution  and  impedance  for  piecewise  - 
straight  wire  loop  antennas.  Results  are  shown  for  planar  examples 
in  this  report;  a  future  report  will  deal  with  the  analysis  of  three- 
dimensional  multiturn  loop  antennas  based  on  the  same  method  reported 
here. 


A  brief  review  of  the  theory  precedes  the  numerical  results. 


II.  THEORY 

Consider  the  polygon  antenna  as  a  piecewise  -  straight  wire  loop 
with  E  and  H  its  radiated  electric  and  magnetic  field  intensities  every¬ 
where  in  space.  The  surface-equivalence  principle  of  Schelkunoff2  is 
used  to  replace  the  antenna  by  equivalent  electric  and  magnetic  surface 
current  densities,  J  and  M,  radiating  in  free  space.  These  quantities 
are  obtained  from  the  relations  J  =  &  X  H  and  M  :  EX  evaluated  at 
the  antenna  surface  with  n  the  unit  normal  outward  from  the  surface. 

The  current  densities  J  and  M  radiate  the  actual  field  intensities  (E,  H) 
in  the  region  exterior  to  the  antenna  and  generate  a  null  field  in  the 
interior  of  the  antenna  itself.  Thus  there  would  exist  no  coupling,  or 
reaction,  between  the  source  distribution  (J,  M)  and  any  source  located 
entirely  within  the  interior  region. 

This  zero- reaction  condition  has  been  utilized  previously*"5  to 
solve  problems  involving  dipole  antennas  and  is  here  extended  to  include 
polygon  loops.  The  loop  is  represented  by  a  set  of  N  tubular  segments, 
each  having  a  piecewise  -  sinusoidal  current  distribution.  If  the  radius 
of  the  wire  is  small  compared  with  the  circumference  of  the  loop,  then 
for  the  purpose  of  calculating  reaction  the  surface  current  density  J 
may  be  replaced  by  a  filamentary  current  I  on  the  surface  of  the  antenna. 

In  this  way  the  set  of  N  tubular  segments  is  replaced  by  a  set  of  fila¬ 
mentary  segments  which  can  be  considered  as  a  collection  of  N  overlapping 
dipoles. 

The  zero-reaction  condition  is  implemented  by  enforcing  reciprocity 
between  the  collection  of  filamentary  surface- current  dipoles  and  a 
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collection  of  N  filamentary  test  dipoles  located  along  the  axis  of  the 
loop.  Figure  1  illustrates  the  situation  for  a  typical  pair  of  elements. 
The  coordinate  Sj  along  the  j-th  surface  dipole  has  unit  vector  8j.  The 
reciprocity  relation  between  the  j-th  true  surface-current  dipole  and 
the  i-th  axial  test  dipole  has  the  following  form: 

(1)  jj  •  Mj  dSj  *  J  Ei  ♦  Jj  dsj 

In  Eq.  (1)  the  currents  are  that  portion  of  the  surface-current 

distribution  associated  with  the  j-th  surface  dipole»  and  (E  H,)  are  the 
free-space  fields  of  the  i-th  test  dipole.  The  surface  integral  on  the 
left  side  of  Eq.  (1)  exists  in  the  vicinity  of  the  j-th  junction  where  Mj 
is  nonzero,  and  the  line  integral  on  the  right  side  is  evaluated  along  the 
j-th  dipole  filament. 


Fig.  1.  Typical  locations  for  i-th  test  dipole  on  axis 
and  j-th  true  surface -current  dipole. 
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After  the  scalar  product  has  been  taken,  the  right  side  of  Eq.  (1) 
normalized  by  the  j-th  terminal  current  Ij  has  the  form  of  mutual 
impedance: 


(2)  zij  =  -  ~  J  Eit(sj)  Ij(sj)  dsj  . 

J 

The  quantity  E^(sj)  represents  the  component  of  the  free-space  electric 
field  radiated  by  the  i-th  test  dipole  that  is  tangential  to  the  j-th  surface 
dipole*  If  a  delta-gap  voltage  generator  (Mj  in  this  case  is  a  ring  of 
magnetic  current)  i:.  assumed  at  the  terminals  of  the  j-th  surface  dipole, 
then  the  left-hand  side  of  Eq.  (1)  reduces  to  zero  except  when  i  and  j 
coincide.  For  this  case,  the  integral  is  equal  to  V^,  the  applied  voltage 
at  that  junction. 

By  enforcing  reciprocity  via  Eqs.  (1)  and  (E)  between  the  i-th  test 
dipole  and  all  N  true  surface-current  dipoles,  a  set  of  N  mutual  imped¬ 
ance  coefficients  is  obtained.  The  coefficients  relate  the  terminal 
voltages  and  currents  as  follows: 

N 

(3)  £  Zij  Ij  =  Vi 


Equation  (3)  represents  one  linear  equation  out  of  a  whole  system  of 
N  equations  that  are  obtained  by  successively  enforcing  reciprocity 
between  the  N  surface  dipoles  and  each  of  the  N  axial  test  dipoles  in 
turn.  This  system  of  equations  is  solved  by  fixing  a  unit  voltage  source 
at  the  feed  terminals  and  short-circuiting  the  other  terminals.  Solution 
of  the  equations  yields  values  for  the  current  samples  at  the  N  junctions, 
thus  determining  the  current  distribution  over  the  entire  perimeter  of 
the  antenna. 

Impedance  information  is  obtained  from  the  current  sample  at 
the  feed  terminals.  Since  a  unit  voltage  source  was  assumed,  this 
quantity  is  itself  the  input  admittance  and  its  reciprocal  is  the  input 
impedance  of  the  antenna. 

In  the  next  section  numerical  data  are  compared  with  experi¬ 
mental  data  and  independent  theoretical  calculations. 


31 


4 


III.  NUMERICAL  RESULTS 


Comparison  with  experiment  and  comparison  with  results  from 
independent  theory  are  two  important  ways  to  demonstrate  the  accuracy 
of  the  multisegment  reaction  technique.  Both  comparisons  will  be  shown 
in  this  section,  along  with  results  for  which  there  are  no  direct  compari¬ 
sons,  but  which  are  nevertheless  interesting,  useful,  and,  it  is  believed, 
accurate . 

Figure  2  compares  measured  and  calculated  admittances  for  half 
a  square  loop  mounted  vertically  over  a  ground  plane.  One  end  of  the 
antenna  is  coaxially  fed  out  of  the  ground  plane.  The  agreement  between 


Fig.  2.  Calculated  and  measured  data  for  the  conductance  (G) 
and  susceptance  (B)  of  a  vertical  half-loop  over  a 
ground  plane. 
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experiment  and  theory  that  is  shown  in  Fig.  2  demonstrates  that  the 
multisegment  reaction  technique  can  be  applied  successfully  to  closed, 
piecewise- straight  wire  antenna  shapes. 

Figures  3-6  are  a  series  of  comparisons  between  the  admittance  of 
a  regular  polygon  and  that  of  a  circular  loop  of  the  same  perimeter.  The 
data  for  the  circular  loop  are  those  of  Storer,  *  who  obtained  his  results 
by  numerically  solving  the  integral  equation  for  a  circular  loop  antenna. 
The  polygon  data  were  calculated  for  symmetrical  side-driven  loops 
using  the  multisegment  reaction  technique. 

It  seems  reasonable  that,  of  all  the  possible  polygonal  shapes 
(for  logical  development,  we  have  chosen  regular  polygons),  the  one 
whose  shape  most  nearly  resembled  a  circular  loop  of  the  same  peri¬ 
meter  would  have  impedance  (or  admittance)  properties  closely  re¬ 
sembling  those  of  the  circular  loop.  Figures  3-6  bear  out  this  reason¬ 
ing.  Notice  that,  as  the  number  of  sides  is  increased,  the  admittance 
curves  more  nearly  resemble  the  circular  loop  data;  for  a  12-sided 
polygon,  the  data  are  virtually  indistinguishable  over  the  entire  range 
of  values.  (Data  calculated  for  a  16-sided  polygon  essentially  lie  on 
that  for  the  circular  loop.)  Notice  also  that  the  electrically  small 
polygon  (P/X  <0.5)  behaves  quite  like  a  circular  loop  for  any  number 
of  sides.  This  bears  out  the  expected  result  -  electrically  small  single¬ 
turn  loops  behave  much  the  same  whether  they  are  circular  or  polygonal 
in  shape. 

Figures  7  and  8  show  that  the  admittance  of  a  planar  polygon  with 
few  sides  is  altered  when  the  feed  point  is  changed  from  a  side  to  a 
corner.  The  curves  in  Figs.  7  and  8  are  for  corner-driven  polygons  of 
three  and  four  sides,  respectively.  These  data  differ  from  the  data  of 
Figs.  3  and  4  most  noticeably  for  higher  values  of  P/  X  .  As  the  number 
of  sides  is  increased,  though,  the  admittances  for  the  corner-driven  and 
the  side-driven  cases  become  nearly  identical.  Thus  the  curves  in  Fig.  5 
and  Fig.  6  apply  to  either  driving-point  situation. 

Figures  9  and  10  demonstrate  that  the  multisegment  reaction 
technique  gives  an  accurate  description  of  the  current  distribution  over 
the  entire  length  of  an  antenna.  The  current  samples  at  junctions  between 
segments  for  a  12-sided,  24-segment  polygon  are  compared  with  data  by 
Storer  for  the  circular  loop  antenna.  The  current  sample  values  com¬ 
pare  well  with  the  circular  loop  data  for  perimeters  of  both  one  and  two 
wavelengths.  The  currents  are  plotted  for  half  of  the  perimeter,  since 
the  current  on  the  other  half  is  symmetric  with  respect  to  the  feeu  ~oint. 
The  sample  data  in  Figs.  9  and  10  represent  either  the  side-driven  or 
the  corner-driven  case,  since  the  difference  in  currents  between  the  two 
cases  at  any  junction  is  less  than  one -tenth  of  one  milliampere. 
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ig.  4.  Conductance  <G)  and  susceptance  (B)  vs.  P/\:  side-driven  four-sided 
polygon  (Richards)  and  circular  loop  (Storer).  (P/<*  =  403.43,  where 


Conductance  (G)  and  susceptance  (B)  vs.  P/\:  side-driven  eight-sided 
polygon  (Richards)  and  circular  loop  (Storer).  (P /a  =  403.43,  where 
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Fig.  7.  Conductance  (G)  and  susceptance  (B)  vs.  P/\:  comer-driven  three- 
sided  polygon  (Richards)  and  circular  loop  (Storer).  (P/o  =  403.43, 
where  a  is  wire  radius.)  * 


Fig.  8.  Conductance  (G)  and  susceptance  (B)  vs  P/\:  corner-driven  four¬ 
sided  polygon  (Richards)  and  circular  loop  (Storer).  (P/ar  =  403.43, 
where  a  is  wire  radius.) 
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Fig.  9.  Current  distribution  on  twolve-sided  polygon  and 
circular  loop:  P/\  =  1.0  and  P/a=  148.41,  where 
a  is  wire  radius,  a  -  real  part  of  current  dis¬ 
tribution;  b  -  imaginary  part  of  current  distri¬ 
bution;  c  -  magnitude  of  current  distribution; 
d  -  phase  of  current  distribution. 
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Fig.  10.  Current  distribution  on  twelve-sided  polygon  and 
circular  loop:  P/\  =  2.0  and  P/o  =  148.41,  where 
(y  ifi  wirp  radius.  a  -  rpal  nart  of  current  d’.Stri 
bution,  b  -  imaginary  part  of  current  distribution; 
c  -  magnitude  of  current  distribution;  d  -  phase  of 
current  dis.ribution. 
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The  above  data  were  calculated  with  a  program  written  in 
Fortran  IV  and  executed  on  an  IBM  7094  digital  computer.  Most  of 
the  computation  is  performed  by  a  part  of  the  program  that  has  been 
written  for  much  more  general  antenna  shapes.  As  an  example  of  the 
time  required,  the  calculation  of  the  input  admittance  and  current  dis¬ 
tribution  for  a  24-segment  polygon  loop  antenna  takes  only  10  seconds. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  multisegment  reaction  technique  has  been  applied  to  calculate 
the  admittance  of  regular  polygon  antennas .  It  has  been  shown  to  give 
good  agreement  both  with  measurements  and  with  independent  theoretical 
calculations.  Moreover,  the  calculated  current  distribution  has  been 
shown  to  be  in  good  agreement  with  that  obtained  from  an  independent 
theory;  hence,  accurate  far-field  patterns  of  polygon  antennas  can  easily 
be  calculated. 

Work  will  continue  on  the  application  of  this  technique  to  three- 
dimensional  polygon  loops  and  multiturn  loops.  It  is  recommended 
that  this  effort  include  a  study  of  how  finite  conductivity  affects  the 
operating  bandwidth  and  radiation  efficiency  of  such  antennas.  The 
use  of  lumped  loading  to  alter  the  current  distribution  alo,.g  an  antenna 
should  be  investigated  for  bandwidth  effects  and  impedance-matching 
possibilities.  Finally,  the  effect  of  a  finite  gap  at  the  feed  point  should 
be  included  to  represent  physical  situations  more  accurately. 
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ABSTRACT 


The  multi  segment  reaction  technique  Is  an  accurate,  efficient 
analytic  tool  for  the  analysis  and  synthesis  of  novel  wire  antenna 
configurations.  The  highlights  in  the  development  of  this  technique 
are  first  presented. 

The  above  technique  is  employed  to  analyze  the  impedance  properties 
of  multi  turn  loops  in  free  space  and  fed  over  a  ground  plane.  A  shift  in 
the  anti  resonant  frequency  is  observed  and  discussed.  The  ability  to  control 
the  anti  resonance  with  tunable  lumped  elements  is  developed.  Finally, 
calculated  and  measured  data  for  square  halo  antennas  of  one  and  two  turns 
are  shown.  Comparison  between  the  two  sets  of  data  is  seen  to  be  good. 


COMPUTER  ANALYSIS  OF  THREE-DIMENSIONAL 
LOOP  ANTENNAS 


I.  INTRODUCTION 

It  is  the  intent  of  this  report  to  relate  results  of  a  computer 
analysis  of  three-dimensional  wire  antennas.  The  particular  shapes 
which  were  analyzed  wore  chosen  because  of  their  similarity  to 
certain  antennas  of  interest  to  the  Ballistic  Research  Laboratory. 

Two  of  these  antennas  are  the  halo- type  antenna  and  the  multi  turn 
loop  antenna.  Previous  efforts  to  improve  the  operation  of  these 
antennas  have  teen  carried  out  on  a  cut-and-try  experimental  basis. 

While  noteworthy  improvements  have  been  made,  it  is  possible  that 
still  better  performance  could  be  achieved  by  implementing  a  com¬ 
puterized  technique  to  sort  through  many  more  alternate  designs  than 
it  would  be  possible  to  test  experimentally. 

A  first  step  toward  a  computerized  approach  of  the  sort  men¬ 
tioned  above  is  the  subject  of  this  report.  The  multi  segment  reaction 
technique  is  a  theoretically  sound  analytic  tool  which  has  been  proven 
on  many  previous  antenna  configurations,  including  straight  and  Vee 
dipoles,  polygon  loons,  and  three-dimensional  antennas  of  various 
shapes.  The  shapes  t.o  be  analyzed  in  this  report  will  further  demon¬ 
strate  the  versatility  of  this  technique.  The  ability  of  the  multi  - 
segment  reaction  t.r '•unique  to  handle  novel  antenna  configurations  and 
to  handle  them  acni.vdy  and  efficiently  is  perhaps  its  greatest  asset 
as  a  tool  for  antenna  design. 

A  very  brief  review  of  the  theory  and  development  of  the  multi¬ 
segment  technique  is  presented  in  the  following  section.  A  later 
section  presents  t.ho  results  of  computer  analysis  along  with  comments 
and  observations. 


A.  The  Multi  segment  Reaction  Technique 

The  theoretical  development  of  this  technique  has  been  presented 
several  times  previous  to  this  report.  For  completeness,  the  important 
points  in  the  development  will  be  repeated  here. 

A  wire  antenna  radiates  the  time-harmonic  fields  (E,H)  in  free 
space.  The  surface  equivalence  theorem  of  Schelkunoff  [l]*is  invoked  to 
replace  the  antonna  by  equivalent  surface  current  densities  (J,M)  which 
radiate  the  fields  (E,H)  exterior  to  the  surface  and  generate  a  null  field 
in  the  interior  region.  Reciprocity  is  enforced  between  the  surface  cur¬ 
rents  and  a  test  source  located  in  the  interior  region,  and  a  zero- 
reaction  [2]  (coupling)  situation  is  seen  to  exist.  The  surface  is  sub¬ 
divided  into  a  collection  of  straight  segments  each  having  an  assumed 
current  distribution  that  is  sinusoidal.  By  considering  the  segments 
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as  a  set  of  overlapping  dipoles  and  by  calculating  the  reaction  between 
each  surface-current  dipole  and  each  member  of  an  equi numerous  set 
of  test  dipoles  all  located  along  the  axis  of  the  wire  antenna,  a 
system  of  linear  equations  is  obtained  which  relates  the  currents  at 
the  junctions  between  each  pair  of  connected  segments  to  the  terminal 
voltages  at  those  junctions.  Assumption  of  an  appropriate  set  of 
voltage  source  conditions  (unit  source  at  driving  terminals,  zero  volts 
at  wire  junctions,  etc.)  results  in  a  solution  of  the  system  of  equations 
for  the  unknown  current  samples.  The  input  impedance  is  obtained  as 
the  reciprocal  of  the  current  sample  at  the  driving  terminal,  and  the 
current  over  the  whole  antenna  is  obtained  in  terms  of  the  sample  point 
currents. 


Multi  turn  Loop  in  Free 


The  multi  turn  loop  antenna  resembles  a  coil  with  several  turns. 
An  individual  turn  may  be  circular,  rectangular,  or  some  other  shape 
and  the  spacing  between  successive  turns  is  generally  small  compared 
to  the  length  of  a  single  turn.  In  free  space  the  multiturn  loop  is 
fed  as  a  balanced  configuration  and  this  presents  a  problem  if  a  coax 
feed  line  is  used.  For  operation  at  a  single  frequency  a  balanced  feed 
may  be  achieved  with  a  tuned  balancing  network  or  a  fixed-length  balun. 
However,  for  operation  over  a  band  of  frequencies  a  wide-band  balun 
of  the  transformer  type  would  be  preferable. 


Figure  1  indicates  the  shape  of  the  loops  with  which  this  section 
is  concerned.  A  two-turn  loop  is  shown,  but  this  analysis  will  cover 
loops  up  to  four  turns  and  spacings  between  turns  of  both  2  cm  and 
4  cm.  For  1 /8-inch  diameter  wire,  this  represents  a  turn  spacing-to- 
wire  radius  ratio  of  about  12.5  and  25,  respectively. 


Figure  2  is  a  plot  of  input  impedance  versus  frequency  (L/x) 
for  a  4-turn  multiturn  loop  with  4  cm  turn  spacing.  The  shape  of  this 
impedance  curve  is  generally  the  same  for  all  antennas  of  this  con¬ 
figuration.  Anti  resonant  peaks  are  found  roughly  at  the  odd-numbered 
half  wavelength  points  (L/x  =  1/2,  3/2,  etc.)  and,  between  peaks,  the 
real  part  of  the  impedance  is  small.  The  number,  relative  width,  and 
location  of  these  peaks  depend  on  both  the  number  of  turns  and  the 
spacing  between  turns.  The  figures  that  follow  demonstrate  how  these 
parameters  affect  the  input  impedance  of  a  multi  turn  loop  antenna. 


One  of  the  most  interesting  phenomena  associated  with  the  per¬ 
formance  of  the  multi  turn  loop  antenna  is  the  influence  of  the  coupling 
between  turns  on  the  frequencies  at  which  anti  resonances  occur.  The 
mechanism  which  induces  this  variation  in  the  location  of  the  peaks  of 
the  impedance  curve  has  been  the  subject  of  discussion  for  some  time. 

It  is  believed  that  a  firmer  understanding  of  this  phenomenon  is  now 
available  through  analysis  of  these  loops  via  the  multisegment  reaction 
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technique.  Figure  3  shows  the  input  impedance  around  the  first  anti¬ 
resonant  point  for  three  square  multi  turn  loops  of  two,  three,  and  four 
turns.  The  turn  spacinq  for  each  loop  is  4  cm.  Note  that  the  value  of 
L/x  at  which  a-. ‘  i resonance  occurs  is  almost  identical  for  all  three 
antennas.  Note  iho  that  the  total  wire  length  is  almost  identical  for 
the  three  antennas.  For  a  given  length  L  of  wire,  the  size  of  each 
individual  turn  decreases  as  antennas  with  more  turns  are  formed  from 
the  wire.  From  this  fiqure  it  is  apparent  that  within  the  range  of 
sizes  investigated,  the  cross-sectional  size  of  the  individual  loop 
(turn)  has  little  if  anything  to  do  with  the  antiresonant  frequency. 

Another  interesting  observation  on  the  data  of  Fig.  3  is  that 
increasing  the  number  of  turns  while  holding  length  constant  does  not 
change  the  anti  resonant  frequency,  given  the  same  turn  spacing.  This 
tends  to  indicate  that  a  balance  is  maintained  between  the  chance  in 
inductance  due  to  a  change  in  the  do forma  tier,  of  the  wire  into  more 
turns  and  the  change  in  the  distributed  capacitance  between  neighboring 
turns  due  to  the  same  deformation.  This  involvement  of  electromagnetic 
coupling  in  the  shifting  of  an  anti  resonant  frequency  (the  thin-wire 
single-turn  loop  of  just  about  any  shape  has  a  natural  antiresonance  very 
near  L/x  =  0.5)  is  admittedly  of  a  complex  nature,  but  these  data  defi¬ 
nitely  indicate  the  presence  of  the  effect  on  the  impedance  of  the  antenna. 

Figure  4  reinforces  the  arguments  made  above.  The  impedance  data 
shown  here  are  concerned  with  two,  three,  and  four  turn  loops  each 
having  a  turn  spacing  of  2  cm.  Again,  essentially  the  same  length 
of  wire  was  used  for  all  three  antennas.  This  length  was  the  same 
as  that  used  for  the  uitennas  in  Fig.  3.  Again,  all  three  antennas 
exhibit  an  anti  resonance  at  nearly  exactly  the  same  value  of  L/x  (same 
frequency).  All  the  statements  made  before  about  the  antennas  of  Fig.  3 
could  be  made  about  ;hese  antennas.  The  only  essential  change  is  that 
the  turn  spacing  has  been  halved  for  the  antennas  of  Fig.  4  from  4  cm  to 
2  cm.  The  inductance  would  not  be  appreciably  affected,  but  the  distributed 
capacitance  between  turns  would  definitely  be  increased.  Notice  that,  for 
the  antennas  in  Fig.  4,  the  antiresonant  frequency  occurs  at  a  value  of 
L/X  near  0.33.  For  the  4  cm-spaced  antennas  of  Fig.  3  the  anti  resonant 
frequency  occurs  at  about  L/x  =  0.38.  Thus  a  reduction  in  the  turn 
spacing  has  a  definite  tendency  to  lower  the  anti  resonant  frequency  of  a 
multi  turn  loop  antenna. 

Figures  5  and  6  present  additional  evidence  that  the  turn 
spacing  really  controls  the  antiresonance  shift.  Figure  5  is  quite  similar 
to  Fig.  3  except  that  the  antennas  in  Fig.  5  all  have  the  same  turn  size 
as  the  3-turn  antenna  of  Fig.  3.  Figure  6  is  like  Fig.  5  except  for 
the  turn  spacing  of  2  cm  instead  of  4  cm.  The  results  are  essentially 
identical  to  those  in  the  corresponding  previous  figures.  The  anti- 
resenant  frequency  for  all  three  antennas  in  Fig.  5  is  very  near  L/x  =  0.38, 
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while  that  for  the  antennas  in  Fig.  6  occurs  near  L/ A  =  0.33.  These 
data  emphasize  the  dependence  of  the  anti  resonant  point  on  the  turn 
spacing,  since  the  turn  size  for  these  latter  cases  is  the  same  for 
all  three  antennas. 

Figure  7  compares  the  impedances  of  t.wo-turn  loops  of  4  cm 
spacing  and  2  cm  spacing.  The  frequency-shifting  effect  of  the  turn 
spacing  on  the  anti  resonant  point  is  less  pronounced  at  the  second 
anti  resonance.  The  situation  is  somewhat  different  in  Fig.  8,  which 
compares  impedance  data  for  four-turn  loops  of  2  cm  and  4  cm  turn 
spacings.  Here  the  loop  with  the  closer  turn  spacing  shows  a  lowered 
second  anti  resonance  as  well  as  the  first.  This  indicates  that  there 
is  still  significant  interaction  between  turns  at  the  higher  frequency. 
This  study,  though,  has  not  dwelt  on  the  impedance  properties  at  higher 
frequencies  for  reasons  which  will  be  given  in  the  following  paragraph. 

The  primary  motivation  for  a  study  of  these  antennas  near  the 
first  anti  resonant  peak  is  that  this  is  the  region  where  the  antenna 
is  ordinarily  designed  to  operate  [3].  It  has  been  seen  that  the  real 
part  of  the  input  impedance  is  low  for  values  of  L/ A  not  near  the  anti- 
resonant  point.  Hence,  if  one  is  to  obtain  a  meaningful  amount  of 
radiated  energy  from  the  antenna,  he  must  operate  the  antenna  near  an 
anti  resonance.  Another  consideration  is  that  of  matching  the  antenna 
impedance  to  the  characteristic  impedance  of  the  feed  system.  Only 
near  anti  resonance  is  there  at  least  some  chance  of  matching  an 
antenna  of  this  configuration  to  a  reasonable  impedance  level,  such 
as  50  ohms  or  100  ohms.  There  is  always  the  possibility  of  tuning  out 
the  reactive  component  of  the  impedance;  however,  more  will  be  said  of 
this  point  in  a  later  section.  One  problem  associated  with  this 
attempt  to  match  impedance  levels  by  operating  near  anti  resonance  is 
the  very  steep  slope  of  the  impedance  curve  in  this  region.  As  one 
adds  turns  (see  Fig.  5)  in  order  to  operate  at  a  lower  actual  fre¬ 
quency,  the  antiresonant  peak  becomes  more  narrow  and  the  slope  is 
even  steeper.  Thus,  the  range  of  frequencies  within  which  a  given 
antenna  performs  satisfactorily  is  quite  narrow,  and  this  range  becomes 
even  more  ne  'ow  as  more  turns  are  added.  This  performance  character¬ 
istic  severe ,j  limits  the  applications  for  which  an  antenna  of  this 
sort  is  useful.  Any  implementation  more  ambitious  than  very  narrow- 
band  operation  would  necessitate  a  rather  elaborate  compensating  system 
due  to  the  rapid  rate  of  change  of  the  antenna  impedance  in  the  operat¬ 
ing  frequency  range.  This  problem  is  not  insurmountable,  however,  and 
material  presented  in  the  next  section  may  be  of  use  in  alleviating 
the  problem. 


C.  Three-Dimensional  Loop  over  a  Ground  Plane 

The  basic  advantage  to  be  gained  by  placing  the  three-dimensional 
or  multi  turn  loop  over  a  ground  plane  is  that  it  can  be  fed  from  a  coaxial 


50 


r*i 


vt 

i 


line  in  an  unbalanced  configuration.  The  presence  of  the  ground  plane  | 

alters  the  behavior  of  the  antenna,  but  if  the  way  in  which  these  j 

alterations  occur  is  known,  then  steps  can  be  taken  to  obtain  l 

satisfactory  performance  from  the  antenna.  Because  the  problem  of  the  1 

balanced  feed  is  overcome,  more  attention  can  be  directed  toward  tun¬ 
ing  the  antenna  itself  to  obtain  improved  radiation  properties.  ; 

Figure  9  is  an  example  of  the  sort  of  antenna  configuration  ; 

with  which  this  section  is  concerned.  The  antenna  in  the  figure  is 
shown  with  2-1/2  turns.  This  discussion  will  include  antennas  from 
1-1/2  to  4-1/2  turns,  with  turn  spacings  of  both  5  cm  and  2  cm. 

Figures  10  and  11  contain  the  information  pertinent  to  this 
discussion.  Figure  10  is  input  impedance  vs  L/x  for  2-1/2  turn  loops 
of  both  5  cm  and  2  cm  turn  spacings.  Figure  11  shows  the  same  data 
for  4-1/2  turn  loops  of  the  same  two  spacings.  (Note:  L  stands  for 
the  length  of  the  antenna  plus  its  image.)  Notice  that,  in  the  case  of 
the  5  cm  turn  spacings,  the  first  antiresonance  actually  occurs  at  a 
value  of  L/x  greater  than  0.5.  This  would  seem  to  indicate  that,  due 
to  the  ground  plane  symmetry,  there  is  a  tendency  to  detune  the  antenna 
so  antiresonance  occurs  at  L/x  above  0.5,  rather  than  below  0.5  as  was 
the  case  for  the  free-space  configuration.  Still,  when  the  turn 
spacing  is  reduced  to  2  cm,  it  can  be  seen  that  the  anti  resonant  point 
drops  below  L/x  =  0.5.  Evidently,  the  proximity  effect  of  nearest- 
neighbor  turns  is  at  least  partially  offset  by  the  presence  of  the 
ground  plane. 

The  use  of  lumped  elements  to  adjust  the  impedance  properties 
of  this  antenna  has  been  probed  only  slightly.  Experimental  work  on  a 
multi  turn  loop  over  a  ground  plane  by  Shreve  [4]  indicates  that  the  im¬ 
pedance  of  a  loop  terminated  in  an  adjustable  capacitance  can  be 
shifted  significantly.  Figure  12  shows  the  net  effect  on  the  impedance 
curve  is  to  shift  the  whole  curve  to  higher  values  of  L/\.  The  amount 
of  shift  is,  of  course,  dependent  on  the  value  of  the  adjustable 
capacitor.  This  illustration  of  tunability  with  regard  to  the  im¬ 
pedance  of  an  antenna  is  valuable  if  one  is  to  consider  operating  the 
antenna  over  any  reasonable  range  of  frequencies. 

Results  of  an  alternate  tuning  technique  are  given  in  Fig.  13. 

Here  the  input  impedance  of  a  1-1/2  turn  loop  is  plotted  against  L/x  for 
several  cases  of  distributed  lumped  capacitances.  In  the  computer 
program  implementing  the  multi  segment  reaction  technique  a  small 
■apacitor  is  inserted  at  every  wire  junction  (the  corners,  in  this 
;ase)  except  the  driving  terminal.  The  impedance  data  in  Fig.  13  show 
a  marked  similarity  to  those  of  Fig.  12.  As  the  value  of  capacitance 
is  increased  (capacitive  reactance  is  decreased),  notice  that  the  im¬ 
pedance  plots  come  closer  to  that  of  the  unloaded  wire  antenna.  Thus 
a  collection  of  adjustable  capacitors  distributed  along  the  antenna  can 
have  much  the  same  effect  on  the  input  impedance  as  a  single  adjustable 


51 


capacitor  at  the  far  end  of  the  loop.  The  advantage  of  the  single¬ 
element  tuning  is  its  simplicity.  The  advantage  Dossessed  by  the 
antenna  with  several  capacitors  is  the  possibility  of  reducing 
voltage  breakdown  [5].  The  capacitors,  in  addition  to  providing  a 
distributed  tuning  network,  may  redistribute  among  themselves  the 
voltage  drop  around  the  loop,  possibly  reducing  a  build-up  of  a 
large  potential  gradient  at  a  single  region  on  the  loop.  This  would 
be  most  useful  if  the  antenna  were  being  operated  at  a  high  level  of 
input  power,  or  if  the  surrounding  medium  contained  ionized  gases. 

Finally,  the  ability  to  match  the  input  impedance  to  a 
desired  value  as  the  frequency  is  changed  is  important.  Figures  12 
and  13  indicate  that,  for  operation  on  the  low-frequency  side  of  the 
anti resonant  peak,  the  reactive  component  of  the  input  impedance  is 
inductive  when  the  real  part  is  around  50  ohms.  Shreve  [6]  has  shown 
that  this  inductive  reactance  may  be  tuned  out  at  the  input  terminals 
with  a  series  capacitor.  Thus,  the  antenna  may  be  continuously  matched 
to  the  feed  system  over  a  significant  frequency  range  through  the  use 
of  adjustable  lumped  elements  with  the  antenna. 

Another  sort  of  loop  antenna  over  a  ground  plane  which  has  re¬ 
ceived  some  attention  is  the  so-called  halo  antenna.  The  models  with 
which  this  study  has  been  concerned  are  only  a  first  step  toward  a 
true  halo  antenna,  but  the  results  which  have  been  obtained  serve  to 
demonstrate  the  accuracy  and  power  of  the  multisegment  reaction  tech¬ 
nique  as  a  tool  both  for  analysis  and  for  synthesis  of  novel  antenna 
shapes. 


Figure  14  shows  a  Smith  chart  impedance  plot  of  calculated  and 
measured  data  for  the  one-turn  square  halo  depicted  in  the  upper 
corner  of  that  figure.  The  total  lengt'  of  the  antenna  plus  its  image 
is  300  cm,  so  the  range  of  values  of  L/a  represented  is  from  L/X  =  2 
to  L/x  =  3.  Notice  that  the  antiresonance  is  delayed,  occurring  at 
a  value  of  L/x  near  2.7,  instead  of  very  close  to  L/x  =  2.5.  This  is 
the  same  effect  noted  previously  when  the  turns  of  the  multi  turn  loop 
were  too  far  apart  for  the  proximity  effect  to  lower  the  anti  resonant 
point  noticeably.  Note  also  the  good  comparison  between  calculated 
and  measured  data. 

Figure  15  illustrates  a  two-turn  halo  antenna  similar  to  the 
one  just  discussed.  Figure  16  is  a  Smith  chart  impedance  plot  for 
the  antenna  of  Fig.  15,  again  comparing  calculated  and  measured  data. 
The  error  between  calculations  and  measurements  at  the  highest 
frequency  shown  is  less  than  2  percent.  The  good  comparison  which 
Fig.  16  shows  between  theory  and  experiment  lends  weight  to  the 
validity  of  the  multisegment  reaction  technique.  Notice  once  again 
that,  although  the  total  wire  length  plus  its  image  equals  the  odd- 
numbered  half-wavelengths  at  frequencies  of  100  MHz  and  300  MHz,  the 
antiresonant  points  are  delayed  in  both  instances.  Again,  this  is 
attributed  to  the  effect  of  the  ground  plane  on  the  impedance  prop¬ 
erties  of  the  antennas  involved  in  this  study. 
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II.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  multi  segment  reaction  technique  has  proven  effective  for 
analyzing  novel  antenna  configurations.  Through  the  use  of  this 
technique  the  input  impedance  of  the  multiturn  loop  has  been  shown 
to  possess  a  predictable  anti  resonance  feature.  The  turn  spacing 
is  a  controlling  factor  in  determining  the  frequency  at  which  this 
anti  resonance  occurs.  The  same  effect  has  been  seen  for  multi  turn 
loops  fed  over  a  ground  plane,  but  the  effect  of  the  ground  plane  it¬ 
self  seems  to  diminish  the  proximity  effect  between  successive  turns. 

The  antenna  has  been  seen  to  be  tunable  by  several  techniques,  two  of 
which  are  the  termination  in  an  adjustable  capacitor  and  the  distributed 
loading  with  several  capacitors.  By  either  means  it  is  possible  to 
adjust  the  frequency  at  which  anti  resonance  occurs,  or  to  realize 
a  given  input  impedance. 

The  results  which  this  report  presents  are  only  the  completion 
of  some  basic  steps  toward  the  realization  of  a  completely  analytic, 
computerized  technique  for  wire  antenna  design.  Yet,  the  accuracy 
with  which  these  calculations  have  been  made,  as  demonstrated  by 
their  comparison  with  experimental  measurements,  shows  that  these  have 
indeed  been  large  steps.  A  greater  understanding  of  the  operation 
of  these  new  antenna  configurations  is  unfolding  via  the  computerized 
analysis.  Further  investigation  of  the  properties  and  possibilities 
of  these  antennas  is  necessary.  More  improvements  ought  to  be  made 
in  the  analytic  technioue  to  allow  more  general  configurations  to  be 
handled.  The  alternatives  of  tunability  through  lumped  loading  should 
be  further  explored. 
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WIRE  RADIUS:  0.1588  cm 


Fig.  1.  Two-turn  square  loop  in  free  space; 
turn  spacing  =  4  cm. 
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Fig.  3(a).  Input  resistance  near  first  anti  resonance: 

Square  loops  of  two,  three,  and  four  turns 
(turn  spacing  =  4  cm,  same  total  wire 
length). 
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L2  *  637 cm 
Lj  s  644cm 
L4  *  643  cm 


WIRE  RADIUS: 
0.1588  cm 


Fig.  4(a) .  Input  resistance  near  first  anti  resonance: 

square  loops  of  two,  three,  and  four  turns 
(turn  spacing  =  2  cm,  same  total  wire 
length). 
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WIRE  RADIUS^ 
0.1588  cm 


TURN  SIZE: 
26cm  x  26cm 


Fig.  6(a).  Input  resistance  near  first  anti  resonance: 

square  loops  of  two,  three,  and  four  turns 
(turn  spacing  =  2  cm,  same  turn  size). 


it  impedances 
spacings  of 


impedances 


FREQUENCY  (MHz) 


Fig.  12(b).  Input  reactance  showing  effect  of  capacitor 
between  multi  turn  loop  and  ground  (Shreve). 


Fig.  13(a).  Variation  of  input  resistance  near  first 

antiresonance  for  several  cases  of  multiple 
reactive  lumped  loads:  one  and  one-half 
turn  vertical  half-loop  over  ground  plane. 


71 


Fig.  13(b).  Variation  of  input  reactance  near  f^rst 

untiresonance  for  several  cases  of  multiple 
reactive  lumped  loads:  one  and  one-half 
turn  vertical  half-loop  over  ground  plane. 


WIRE  RADIUS  :  o  =  0.1588  cm 


Fig.  14.  Calculated  and  measured  input  impedance  of 
one- turn  halo  antenna. 
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WIRE  MATERIAL  :  BRASS 
WIRE  RADIUS  :  a  =  0.I58B  cm 
WIRE  LENGTH  I  L  =  75  cm 


i  it  »wv«  ~».v»  m 


,  <‘'-fft*n*‘r**W •rMtSX&X# 


300  MHz  | 


200  MHz' 


■  *— •  CALCULATED 
o  o  MEASURED 


Z0*  50  ohms 


Fig.  16.  Calculated  and  measured  input  impedance  of  two-turn 
halo- type  antenna  over  ground  plane. 
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ABSTRACT 


The  Induced  emf  formulation  is  employed  to  develop  a  closed  form 
expression  for  the  mutual  impedance  between  coplanar-skew  dipoles. 
Numerical  results  are  presented  in  graphical  form. 


MUTUAL  IMPEDANCE  BETWEEN  COPLANAR-SKEW  DIPOLES 


I .  INTRODUCTION 
1  * 

H.E.  King  has  determined  the  mutual  impedance  between  parallel 
dipoles,  and  Lewin?  has  analyzed  the  coplanar-skew  dipoles.  Lewin's 
analysis,  however,  is  restricted  to  half-wave  center-fed  dipoles.  In 
this  paper  we  derive  an  expression  for  the  mutual  impedance  between 
coplanar-skew  dipoles  with  arbitrary  lengths  and  terminal  positions. 

This  expression  *:s  relatively  simple  and  convenient  for  computer 
programming.  Future  reports  will  use  this  result  to  develop  a  multi  - 
segment  reaction  formulation  for  multi-turn  loops  and  halo  antennas. 

In  the  following  section  we  present  a  new  expression  for  the  near¬ 
zone  field  of  a  linear  dipole,  in  a  form  most  convenient  for  the  mutual- 
impedance  analysis. 


II.  THE  NEAR-ZONE  FIELD  OF  THE  LINEAR  DIPOLE 

As  shown  in  Fig.  1,  we  consider  a  dipole  located  on  the  z  axis 
and  let  and  Z3  denote  the  endpoints  and  Z2  the  terminals.  In  the 
induced  emf  method,  the  dipole  current  is  given  by 

sin  k(z-z,) 

(1>  I(z)  =  2  'l  sin  kc l  for  21  <  2  <  z2 

sin  k(z--z) 

(2>  I<2>  *  2  !1  ~  sin'  kc'2  ■  for  z2  <  2  <  z3 


The  time. dependence  e^“2  is  suppressed,  cj  and  C£  denote  the  dipole 
arm  lengths,  z  is  a  unit  vector,  I i  is  the  terminal  current,  and  k  =  2i/\. 
The  field  generated  in  free  space  by  this  dipole  (or  sinusoidal  line  source) 
is  determined  most  readily  from  the  expressions  of  Schelkunoff .3  The  cy¬ 
lindrical  components  of  the  electric  field  intensity  are 


(3)  E_ 


r  -jkR, 

j  301  j  (z-z^e 

p  Rj  sin  kCj 


-jkR2  "JkR3 

(z-z2)e  sin  kc  (z-z3)e 

R2  sin  kc^  sin  kc^  +  R3  sin  kc2 


(4)  Ez  =  j30I x 


-jkR, 

e  ,  .  + 

sin  kCj 


-  jkRp 

e  sin  kc 
r’2  sin  kc^  sin  kc2 


-jkR~ 

e _ 

R3  sin  kc2 


* 

References  are  listed  on  p.  85. 
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where  c  is  the  dipole  length  and  R.,  R0  and  R,  are  the  distances  defined 
In  Fig.  1.  i  Z  i 


Fig.  1.  A  linear  dipole  and  the  coordinate  system. 


The  radial  component  of  the  field  is  obtained  from  a  linear 
combination  of  E  and  E  : 

p  z 

-J301 ,  3  ■3kRm 

(5)  E„  «  — —  f  c_  z_  s-j - 

r  r  m  ra  "m 


where 
(6)  C,  =  r? 


1 


1  sin  kc 


1 
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^  c3  B  Sinkc2 

Equation  (5)  Is  believed  to  be  a  new  and  useful  form  for  the  rigorous 
near- zone  field  of  the  linear  dipole. 


III.  THE  MUTUAL  IMPEDANCE 

For  the  coplanar-skew  dipoles  shown  in  Fig.  2,  the  mutual  impedance 
Is 

r3 

<9>  Z12  *  '  TTTT  [  J»<r>  '  -l(r)dr 

rl 


j 
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The  current  on  dipole  2  Is 


.  I?  sin  k(r-rj 

(10)  i«(r)  ■  r  —  sin  kd,  '  for  rl  <  r  ‘  r2 


(11)  yr)  -  r 


I2  sin  k(r3-r) 


for  r2  <  r  <  r3 


From  Eqs.  (5)  and  (9), 

130  3  i""3  dr 

(12)z^fJi2^J  ^7— 


These  Integrals  are  given  In  terms  of  the  sine-  and  cosine-integrals 
as  follows 

tr2  oijr  -jzm  jzm 

(13)  zm  - —  =  [e  m  E(R  -  z  +  r)  -e  m  E(R  +zm+r)  ] 

m  r  Rm  *■  mm  mm, 

m 


where 

(14)  E(x)  =  Cl ( | x 1 )  -  j  Si (x) . 


From  Eqs.  (12)  and  (13),  the  mutual  impedance  is 


33  11  jk(pz  +qr  ) 

(!5)  Z12  -  -!5  l  l  CD  l  l  pqe  m  " 

m=l  n=l  m  n  p=-l  q=-l 


E(kRmn  +  kPVk^n} 


where  p  and  q  assume  only  the  values  ±1  and  Rmn  is  the  distance  from 
point  zm  on  dipole  1  to  point  rR  on  dipole  2: 

<16)  Rmn  =  *  Tl  '22m  rn  cos  *  ' 


The  coefficients  D„  have  the  same  form  as  the  Cm: 

n  m 


(17)  Dj  =  jrl 


.82 


(18) 


D  *•  -  sin  kd 
u2  sin  kdj  sin  kd2 


<19>  °3  ■  dhag 


The  coordinates  zm  and  rn  are  measured  from  the  coordinate  origin 
at  the  intersection  of  the  axes  of  the  coplanar  skew  dipoles.  For  center 
fed  half-wave  dipoles,  Eq„  (15)  reduces  to  Eq.  (15)  of  Lewin's  paper. 2 

For  parallel  dipoles,  zm  and  rn  go  to  infinity  and  it  is  not  diffi¬ 
cult  to  show  that  Eq.  (15)  reduces  to  the  expression  published  by 
H.E.  King.1 


IV.  NUMERICAL  RESULTS 

Figure  3  illustrates  the  mutual  impedance  between  half-wave  center- 
fed  dipoles  as  a  function  of  the  angle  between  their  axes.  Using 
Eq.  (15),  the  calculations  require  0.15  seconds  on  an  IBM  7094  computer 
for  each  value  of  >i>. 

Although  the  numerical  approach  of  Baker  and  LaGrone^  is  efficient 
when  the  dipoles  are  adequately  separated,  Eq.  (15)  is  preferable  for 
closely  spaced  dipoles. 


V.  CONCLUSIONS 

A  convenient  new  form  is  presented  for  the  rigorous  near-zone 
field  of  the  linear  dipole  with  sinusoidal  current  distribution.  This 
is  employed  to  derive  a  closed-form  expression  for  the  mutual  impedance 
between  coplanar-skew  dipoles.  Future  reports  will  use  these  results 
to  develop  a  multi -spvsf'ot  reaction  formulation  for  halo  antennas  and 
multi -turn  loons. 
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Fig,  3.  Mut  i  impiuance  between  coplanar 
skt  dipoles. 
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Preceding  page  blank 


ABSTRACT 


The  reaction  concept  is  applied  to  determine  the  self-  and 
mutual- admittances  of  coupled  linear  antennas  with  arbitrary  positions 
and  angular  orientations.  Numerical  results  are  presented  in  graphical 
form  for  parallel,  coplanar-  skew,  and  nonplanar-skew  cylindrical-wire 
antennas . 
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COUPLED  LINEAR  ANTENNAS  WITH  SKEW  ORIENTATION 


I.  INTRODUCTION 

i  * 

Chang  and  King  have  solved  the  coupled-dipole  problem  for  the 
special  case  where  the  antennas  are  parallel.  In  this  report  we  con¬ 
sider  the  skew  array  of  linear  antennas.  Future  reports  will  extend 
the  techniques  to  various  three-dimensional  antennas  such  as  the 
multi-turn  loop  and  the  halo. 

The  theory  is  developed  in  the  next  section,  followed  by  numerical 
results  for  some  parallel  and  skew  configurations. 


II.  THEORY 


Consider  the  coupled  linear  antennas  shown  in  Fig.  1.  These 
antennas  are  assumed  to  be  perfectly  conducting,  but  they  may  be 
parallel,  coplanar  skew,  or  non-planar  skew.  The  wire  is  divided 
into  segments,  and  the  complex  constants  Ij,  I2 ,  .  .  Ipj  denote  the 
current  function  sampled  at  the  junctions  of  these  segments.  If  the 
radius  of  these  cylindrical-wire  antennas  is  small,  we  assume  the 
current  vanishes  at  the  ends.  Furthermore,  the  current  density  is 
assumed  uniform  around  the  circumference  of  each  wire;  this  is 
reasonable  if  the  spacing  between  them  is  not  too  small. 


if  antenna  1  is  divided  into  four  equal  segments  as  illustrated  in 
Fig.  1,  a  convenient  expansion  for  the  current  is 


(1) 


Pn(s) 


sinkf 


where  l  is  the  segment  length,  k  =  2t t/\,  and  s  is  the  coordinate  distance 
measured  along  the  axis  of  antenna  1.  The  pulse  function  Pn(s)  has  unit 
value  for  sn-i  <  s  <  sn+i  and  vanishes  elsewhere.  The  unit  vector  s 
defines  the  reference  direction  for  positive  current  flow.  The  time 
dependence  eJwt  is  assumed  and  suppressed. 


The  current  on  antenna  2  is  given  in  a  similar  manner  by 


(2) 


A 

I (t)  =  t 


In  PnW 


sin k(P  -  |t  -  tn| ) 
sinkf 


References  are  listed  on  p.  IOC. 
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Fig.  1.  The  coupled  antennas  are  divided  into  segments 
for  the  piocewise-sinusoidal  current  expansion. 


The  segment  length  £'  on  antenna  2  may  differ  from  that  on  antenna  1. 

Of  course,  Eqs.  (1)  and  (2)  are  not  rigorous  expressions  for  the 
current  functions.  However,  they  represent  current  distributions 
which  are  continuous  and  vanish  at  the  ends  of  each  antenna.  Further¬ 
more,  they  allow  for  slope  discontinuities  across  any  generators  or 
lumped  loads  located  at  the  current  sampling  points.  A  tradeoff  between 
accuracy  and  computational  expense  can  be  achieved  simply  by  dividing 
each  wire  into  a  suitable  number  of  segments.  Our  objective  is  to  deter¬ 
mine  the  expansion  coefficients  or  the  current  samples  In  and  thereby 
the  self-  and  mutual- admittances  of  the  coupled  antennas. 

It  is  helpful  to  apply  the  surface -equivaler.ee  principle  of 
Schelkunoff2  to  the  coupled  antenna  problem.  Let  S  denote  the  circular- 
cylindrical  surfaces  ot  the  antennas.  To  make  S  a  closed  surface,  we 
include  the  flat  ends  of  each  antenna.  Let  E  and  H  represent  the  electric 


90 


and  magnetic  field  intensities  generated  by  the  antenna  array.  Then 
the  equivalent  electric  current  density  is  J  =  n  X  H  and  the  equivalent 
magnetic  current  density  is  M  =  E  X  n.  These  currents,  located  on  S 
and  radiating  in  free  space,  generate  the  true  field  (E,  H)  in  the 
exterior  region  and  a  null  field  in  the  interior.  The  outward  unit  normal 
vector  on  S  is  denoted  by  ft. 

Since  J  and  M  generate  a  null  field  in  the  interior  region  of  the 
antennas,  they  have  no  coupling  or  reaction  with  any  source  located 
in  the  interior.  The  reaction  with  an  electric  line  source  I'(s)  or  I'(t) 
on  the  axis  of  antenna  1  or  antenna  2  is  given  by 


(3) 


-  M  •  H')  dS  =  0 


where  (E'.H1)  denotes  the  free-space  field  generated  by  the  interior 
line  source.  Equation  (3)  follows  from  the  reciprocity  theorem  and 
the  fact  that  J  and  M  generate  a  null  field  in  the  interior.  The  current 
density  J  is  related  to  the  current  I  as  follows 


I 


where  a  denotes  the  wire  radius. 


Let  us  enforce  Eq.  (3)  with  a  sinusoidal  electric  dipole  located 
on  the  axis  of  antenna  1  or  antenna  2.  We  assume  there  are  no  lumped 
loads  and  there  is  only  one  generator:  a  one- volt  generator  at  the 
center  of  antenna  1 .  In  the  first  reaction  test,  the  electric  dipole 
extends  from  sQ  to  on  the  axis  of  antenna  1 .  In  the  second  reaction 
test,  the  dipole  extends  from  sj  to  S3.  In  reaction  tests  4,  5,  and  6,  the 
test  dipole  is  located  <  •  various  positions  on  the  axis,  of  antenna  2,  and 
in  reaction  test  m  tin  Id  on  the  test  dipole  is  denoted  (Em ',  Hm‘) .  If 
the  test  dipole  has  un  1  .  <  rminal  current,  Eq.  (3)  reduces  to 


5. 


Pm'  ds 


-f 


l(t)  •  Ejn'  dt 


=  V 


m 


In  Eq.  (f>),  the  field  Ejr)'  is  evaluated  on  the  surface  of  antenna  1  or  2 . 
Vm  has  unit  value  if  m  =  2  and  is  zero  otherwise.  As  stated  earlier, 
it  is  assumed  that  the  wire  radius  "a"  is  small.  Furthermore,  we 
are  using  a  "delta-gap"  model  for  the  voltage  generator  although  it  is 
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not  difficult  to  develop  the  theory  for  a  more  realistic  generator  with 
finite  extent* 


From  Eqs.  (1),  (2)  and  (5)  we  find  that 

N 

(8)  ^  In  Zmn  ■  Vm 

n=l 

where  Zmn  is  the  mutual  impedance  between  the  axi.ul  test  dipole  m 
and  the  surface  dipol  ■  ,  node  n.  Specifically, 


rrin+i 

(7) 

Zmn  ~  - 

*>11-1 

and 

pfn+l 

(8) 

Zmn  =  ■ 

[ 

tn-i 

A  sink(f  -  I  s>  -  sn  ) 

s  •  Ern  (s)  - - — -  ds  if  1  <  n  <  3 

— ‘  smk£ 


t  •  E'(t) 


sink(£'  -  t  -  tn  ) 


sin  kf 


dt  if  4  <  n  <  N  . 


Dipoles  m  and  n  are  filamentary  dipoles  with  sinusoidal  current  dis¬ 
tributions.  Thus,  the  elements  in  our  mutual-impedance  matrix  are 
given  by  classical  induced  emf  theory.  If  dipoles  m  and  n  are  located 
on  the  same  antenna,  they  form  essentially  a  collinear  array,  but  with 
a  slight  displacement  equal  to  the  wire  radius  a.  Then  Z^n  can 
expressed  in  terms  of  sine  integrals  and  cosine  integrals.  If  dipoles 
m  and  n  are  located  on  opposite  antennas,  the  numerical  technique  of 
Baker  and  LaGrone4  is  applicable.  Alternatively,  if  antennas  1  and  2 
are  coplanar-skew,  the  closed-form  expression  of  Richmond  and 
Geary5  may  be  employed. 

Letting  N  represent  the  total  number  of  unknown  current  samples 
on  our  two- antenna  array,  we  enforce  N  reaction  tests  by  letting  m  =  1, 
2,  . .  N  in  Eq.  (6).  This  yields  a  system  of  simultaneous  linear 
equations.  This  system  is  solved  numerically  to  determine  the  current 
samples  In.  The  impedance  matrix  ZrnTl  is,  of  course,  symmetric. 

The  self  admittance  is  Yn  =  I*  and  the  mutual  admittance  is  Yj2  =  I5  . 

If  each  antenna  is  divided  into  just  two  segments,  our  solution 
reduces  to  the  classical  induced  emf  theory  and  the  accuracy  is  limited. 
As  each  antenna  is  subdivided  into  more  segments,  our  solution  con¬ 
verges  rapidly  and  shows  excellent  agreement  with  experimental 
measurements.  (An  exception  is  the  self-susceptance  Bn  which  has 


poor  convergence  with  the  delta- gap  model.  This  problem  disappears, 
however,  when  the  finite  extent  of  the  generator  gap  is  taken  into 
account.) 


III.  NUMERICAL  RESULTS 

Figures  2,  3  and  4  illustrate  our  numerical  results  for  the  self- 
and  mutual- admittance  versus  the  displacement  d  for  an  array  of 
parallel  nonstaggered  antennas.  Each  antenna  is  divided  into  ten  seg¬ 
ments  for  the  reaction  solution.  It  may  be  noted  that  the  results  show 
excellent  agreement  with  the  data  published  by  Chang  and  King.1 

Figures  5  and  6  show  numerical  results  for  the  mutual  admittance 
between  ccplanar-skew  antennas  versus  the  angle  4*  between  their  axes. 
When  ip  =  0  the  antennas  are  parallel  and  nonstaggered,  and  when 
ij;  =  180*  they  are  collinear.  In  these  limiting  cases  our  results  agree 
closely  with  those  of  Chang  and  King.1 

Figure  7  shows  the  mutual  admittance  between  nonplanar-skew 
antennas  versus  the  angle  <£•  When  <{>  -  0  the  antennas  are  parallel 
and  nonstaggered.  When  4>  s  90®  they  are  orthogonal  and  the  mutual 
admittance  vanishes. 

The  computer  program  developed  for  this  study  handles  linear 
antennas  wiv.h  arbitrary  lengths  Li  and  L2,  displacement  d,  stagger  s 
and  angles  i|j  and  The  output  data  include  the  current  distributions, 
radiation  patterns  and  radiation  efficiency.  Using  an  IBM  7094  com¬ 
puter.  the  computation  time  is  four  seconds  for  each  coupled  antenna 
problem . 


IV.  CONCLUSIONS 

A  multi- segment  reaction  technique  is  developed  and  applied  to 
coupled  linear  antennas.  Numerical  results  are  presented  in  graphical 
form  for  parallel,  coplanar-skew  and  nonplanar-skew  cylindrical-wire 
antennas.  When  the  antennas  are  parallel  or  collinear,  the  calculated 
results  show  excellent  agreement  with  previously  published  data. 

Future  reports  will  extend  these  techniques  to  various  three- 
dimensional  antennas  such  as  the  multi-turn  loop  and  the  halo. 
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Fig.  2.  Self  admittance  vs.  spacing  for  parallel 
half-wave  antennas. 
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Fig.  3.  Mutual  admittance  vs.  spacing  for  parallel 
half-wave  antennas . 


Fig.  4.  Mutual  admittance  vs.  spacing  for  parallel 
one -wave  antennas. 
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Fig.  5.  Mutual  admittance  vs.  rotation  angle  for 
half-wave  coplanar-skew  antennas. 
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.  Mutual  admittance  vs.  rotation  angle  for 
one-wave  coplanar-skew  antennas.  ’ 
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Fig.  7.  Mutual  admittance  vs.  rotation  angle  for 
half-wave  nonplanar-skew  antennas. 
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ABSTRACT 


The  reaction  concept  is  employed  to  formulate  n  solution  for  thin- 
wire  antennas  with  arbitrary  shape.  Based  on  this  theory,  a  digital 
computer  program  has  been  developed  to  provide  a  complete  analysis  of 
such  antennas.  Calculated  and  measured  impedances  are  included  for 
several  antennas  to  demonstrate  the  accuracy  and  versatility  of  the 
formulation  and  t  u  program. 

This  report  presents  the  computer  program  and  the  details  of  the 
solution  for  the  current  distribution,  impedance,  far-zone  patterns, 
radiation  efficiency  and  scattering.  The  analysis  of  wires  with 
multiple  excitat ion  and  lumped  loads  is  considered. 

The  techniques  can  be  extended  to  a  more  advanced  problem  of 
considerable  current  interest:  A  wire  antenna  exciting  a  circumfer¬ 
ential  slot  on  a  rocket  or  projectile. 
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I.  INTRODUCTION 


The  upper  atmosphere  research  program  of  the  Signature  and 
Propagation  Division  of  the  Ballistic  Research  Laboratory  presents 
severe  design  problems  in  the  field  oi  rocket  mounted  antennas. 

Space,  weight,  and  .-onf  iguration  limitations  Impose  severe  restric¬ 
tions  on  design.  M,:h  radiating  efficiencies  and  uniform  field 
strength  patterns  are  required  to  make  accurate  measurements  of  the 
parameters  under  investigation  as  they  are  affected  by  both  natural 
and  artificial  dinurbances.  Because  of  the  small  physical  size  of 
the  racket  vehicle  the  antennas  are  necessarily  short  compared  to 
those  normally  used  for  the  frequencies  when  r.pace  is  not  a  problem. 

To  prevent  high  voltages  across  the  antennas  and  the  resultant  voltage 
breakdown  and  low  efficiency  the  standing  wave  ratio  at  the  antenna 
feed  point  must  be  kept  low.  The  Polar  Cap  Absorption  experiment 
presents  a  typical  problem.  Six  antennas  operating  at  frequencies  of 
9,  18,  36,  72,  144  and  576  MHz  must,  in  effect,  occupy  the  same 
limited  space  on  a  nine  inch  diameter  rocket.  At  the  9  MHz  frequency, 
for  example,  the  maximum  dimension  of  the  effective  radiating  element 
can  be  only  0.068  wavelengths.  This  imposes  very  severe  design  and  oper¬ 
ational  problems. 

Novel  antenna  configurations  which  seem  to  provide  greatly 
improved  operational  characteristics  have  recently  been  devised.  How¬ 
ever,  the  progress  to  date  has  been  primarily  by  the  application  of 
intuitive  "cut  and  try"  techniques  rather  than  the  mere  reliable 
methods  of  analysis  based  on  an  understanding  of  the  radiating 
mechanism  and  the  relations  between  radiation  pattern,  polarization, 
impedance,  bandwidth,  and  other  parameters  which  can  be  derived  from 
a  study  of  computerized  mathematical  models  of  these  novel  antennas. 

To  analyze  and  optimize  various  antennas,  The  ElectroScienca 
Laboratory  has  developed  a  computer  program  for  thin-wire  antennas 
mounted  on  a  ground  plane  or  in  free  space.  This  program  will  handle 
a  wide  variety  of  artenna  configurations  including  the  linear  dipole, 
V-dipole,  zig-zag  and  three-dimensional  antennas  such  as  the  halo  and 
multiturn  loop.  The  output  data  include  the  current  distribution, 
impedance,  radiation  efficiency,  directivity  and  far-zone  field 
patterns , 

The  theory  is  developed  in  the  following  sections,  followed  by 
a  comparison  of  calculated  and  measured  results  for  several  antenna 
configurations.  Appendices  present  the  computer  program  and  an 
analysis  of  the  scattering  problem. 


II.  NEAR-ZONE  FIELD  OF  SINUSOIDAL  LINE  SOURCE 

Our  formulation  for  three-dimensional  wire  antennas  uses  a 


103 


piecewise-sinusoiunl  expansion  for  the  current  distribution.  Consider 
a  typical  straight  segment  of  the  antenna  to  be  located  on  the  z  axis 
with  endpoints  at  z,  and  z2  as  shown  in  Fig.  1.  The  assumed  current 
distribution  on  thi3  segment  is 

(1)  I(z)  *  Ae^z  +  Be“^kz 


•Fig.  1.  An  electric  line  source  on  the  z  axis  and 
the  observation  point  at  (p,z). 

where  A  and  B  are  complex  constants,  k  •  2n/A  and  the  time  dependence 
e-*5  is  understood.  The  field  of  this  sinusoidal  line  source  is 
given  by  Schelkunoffl  as  follows 


■t/t 


(3)'  ^ujwcp  Ep  -  U2  CO802  -Jkl2]  e”^kR2  -  [ij  cosgj  -Jklj  ]  e“^Rl 

vhere  Ii  and  I2  denote  I(z.)  and  I(zj,) »  and  ij  and  I2  denote  the 
derivative  of  I(z)  at  zjana  z2.  The  remaining  quantities  are  defined 
in  Fig.  1. 

The  constants  A  and  B  in  Eq.  1  can  be  eliminated  to  express  the 
current  distribution  in  terms  of  the  endpoint  currents  Ijand  I2  as 
follows : 


(4)  I(z)  «  [Ij  sin  k(z2-z)  +  I2  sin  k(z-zj)  ]  /sin  kd 

where  d  is  the  length  of  the  segment.  From  Eq.  4, 

(5)  Ij  «  k(I2-Ij  cos  kd)/sin  kd 

(6)  I2  *  k(I2  cos  kd  -  Ij  )/sin  kd 

from  Eqs.  2-6,  the  electric  field  intensity  is 

(7)  E  -  ^  [(i,  cos  kd-I2)  e~JkRI  +  (I2  cos  kd-Ij)  e"JkR21 

sin  kd  Ri  R* 

<8>  ^  Y*sln“kd  cos  kdXz-ZjJe^1  -  j  Ij  sin  kd  e  ^ 

p  Ki 

+(Ij-I2  cos  kd)(z-z2)  e  ^kRz  +  j  I2  sin  kd  e  ^RRz  ^ 

R2 

Nov  consider  the  radial  field  Er  as  shown  in  Fig.  1.  From 
Eqs.  7  and  8 

(9)  Br  "rasin'  kd  ^*1  cos  kd_I2>  zi  “  J  Ii  sin  kd  e_JkRl 

R1 

+  (I2  cos  kd  -Ij)  z2  e  ^RR2  +  j  I2  sin  kd  e  ^  2J 

R2 
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Eqs.  7-9  are  useful  in  calculating  the  mutual  Impedances  among 
the  various  portions  of  a  wire  antenna. 

The  next  section  shows  how  the  three-dimensional  wire  antenna 
problem  is  reduced  to  a  simpler  one:  the  mutual  impedance  between 
elementary  (two-segment)  V-dipoles. 


III.  THEORY  OF  THE  THREE-DIMENSIONAL  WIRE  \:."lhNNA 

As  indicated  in  Fig.  2,  the  wire  antenna  is  assumed  to  have 
straight  segments  with  lengths  d^  d2,  ••♦dfJ+1  .  Starting  from  one 
end,  the  distance  along  the  antenna  ir.  measured  by  the  coordinate  H. 
The  current  dist; .  .mtion  1(H)  is  unknown.  However,  a  series  of 


/ 


Fig.  2.  A  three-dimensional  thin-wire  antenna  with 
straight  segments. 
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complex  constants  Ij,  I2,  are  defined  to  denote  consecutive 

samples  of  the  current  at  the  points  £1,  lz,  **’!«•  f°r  example, 

II  ■  Klj)  and  I2  *  I(£2)  38  shown  in  Fig.  2.  The  current  is  assumed 
tb  vanish  at  both  ends  of  the  antenna. 

For  accurate  results,  the  segment  length  should  be  about  X/4 
or  smaller.  If  any  straight  segment  of  the  antenna  has  a  length 
greater  than  this.  It  is  subdivided. 

Using  the  piecewise-sinusoidal  approximation,  the  current 
distribution  on  the  antenna  is  given  by 

(10)  i(0  -  l  i  UnPnU>  8in  k<a~V.i) 

n-1  n  sin  kd 

n 


+  sn+l  P„+l<*> 


sin  k(£n+1>jl) 

sin  kd  , , 
n+1 


] 


where  IR  *  I(S.n)  a°d  dR  is  the  length  of  segment  n.  The  unit  vector 
8n  is  parallel  with  segment  n  and  is  directed  from  the  point 
toward  the  point  fcn.  Pn(D  is  a  pulse  function  which  has  unit 
value  on  segment,  n  and  vanishes  elsewhere.  Eqs.  10  represents  a 
current  distribution  which  is  continuous  and  vanishes  at  the  ends  of 
the  antenna.  Furthermore,  this  expression  allows  for  a  slope  dis¬ 
continuity  in  the  current  function  at  each  generator.  Our  objective 
is  to  determine  the  complex  constants  In. 

Let  S  denote  the  composite  circular-cylindrical  surface  of 
the  antenna.  To  make  S  a  closed  surface,  we  include  the  flat  ends. 

As  indicated  in  Fig.  Z,  consider  the  antenna  to  be  excited  by  a 
source  with  electric  and  magnetic  current  densities  and  Mj:  located 
in  an  exterior  region  V.  The  incident  field  (E^H^)  is  generated  by 
(Jt»Hi)  when  radiating  in  free  space.  The  total  field  (E,,H)  is 
generated  by  when  radiating  in  the  presence  ..of  the  wire.  The 

equivalent  surface-current  densities  on  S  are  ■  n*  H  and  «=  E*n, 
where  n  is  the  outward  unit  normal  vector  on  S.  Consider  the  field 
generated  in  free  space  when  the  primary  source  (J^.M^)  and  the 
surface  currents  (Jg,^)  are  both  turned  on.  According  to  the 
surface  equivalence  principle  of  Schelkunoff,  the  total  field 
(E,H)  will  exist  in  the  exterior  region  and  a  null  field  in  the  interior 
region  of  S. 

Since  these  sources  generated  a  null  field  in  the  interior 
region,  they  have  no  coupling  or  reaction  with  any  source  located 
in  the  interior.  The  reaction  with  an  interior  source  is  given  by 
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(11) 


#8  Qs-I-V^ds +  /// y  Orl'-Ht-n’)  dv  -  0 


where  (JE' ,H')  denotes  the  free-space  field  of  the  interior  test 
sour.ce.  Eq.  11  follows  from  the  reciprocity  theorem  of  Carson. ^ 

The  volume  integral  in  Eq.  11  may  reduce  to  a  surface  or  line 
Integral  depending  on  the  nature  of  the  source 

To  simplify  the  first  integral  in  Eq.  11,  we  assume  the  wire 
antenna  has  perfect  conductivity  so  that  Mg  vanishes.  Furthermore, 
we  neglect  the  integration  over  the  flat  ends  of  the  antenna.  To 
simplify  the  last  integral,  let  *  0  and  assume  the  antenna  is 
driven  by  a  series  of  voltage  generators  Vj,V2,***V^  located  at 
£l,  ®>2 *  ***  (A  generator  voltage  is  considered  positive  if  it 

tends  to  force  a  current  flow  in  the  direction  of  increasing  £.) 

For  an  interior  test  source  we  choose  an  electric  line  source  I'U) 
located  on  the  antenna  axis  and  radiating  in  free  space.  Now  Eq.  11 
reduces  to 

<i2)  -//,  a ,•£’■>*  -  ?  \  i'<y 

n»l 

In  deriving  Eq.  12,  the  wire  radius  "a"  is  assumed  to  be  small  and  a 
"delta-gap"  model  is  employed  for  the  voltage  generators.  Thus, 
each  generator  is  represented  by  a  magnetic  current  loop  encircling 
the  wire.  Furthermore,  we  may  use  Ampere's  law: 

,2tt  ,  , 

(13)  a  fQ  (a,£)  d$  =  i  (Z) 

Eq.  13  neglects  the  displacement  current.  (Although  it  is  easy  to 
treat  the  gap  problem  in  a  more  satisfactory  manner,  this  subject  is 
reserved  for  a  future  report.) 

The  current  density  Jg  is  related  to  rhe  current  1^  as  follows 


(14)  4  “  KO 
7ita~ 


From  Eqs.  12  and  14, 

L  N 

(15)  -/  I(£)-E’d£  -  l  Vn  l'(ln) 

0  n*l 

where  L  is  the  overall  length  of  the  antenna. 
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We  shall  enforce  the  reaction  test  (Eq.  15)  with  several 
different  test  sources.  In  test  m,  the  interior  test  source  is 
an  electric  dipole  with  a  sinusoidal  current  distribution  given  by 

(16)  i;  U)  -  amPm(£)  sin  k  (l-lari)  +  s^  P^U)  sin  kQ^-R) 


sin  kd; 


m 


sin  kd. 


rofl 


This  test  source  is  generally  a  V  dipole  (rather  than  a  linear  dipole), 
and  it  nay  be  noted  in  Eq.  16  that  it  has  unit  terminal  current. 

The  test  dipole  is  located  on  the  antenna  axis  and  extends  from  , 

to  £nrt-i*  When  radiating  in  free  space,  it  generates  the  field  E^. 
From  Eqs.  15  and  16, 


(17)  -  / I<*)‘]C  d£  -  Vm 


since  I'(£n)  a  0  and  I'(£m)  “  1.  From  Eqs.  10  and  17, 


<18>  I  h  Zmn  •  V, 

rrl 


m 


where  Z  is  the  mutual  impedance  between  the  axial  test  dipole  m  and 
the  surface  dipole  mode  n.  Specifically, 


<»>  zm„  -  -  / 


tn-l 


sin  k(£-;.n_i )  d£ 
sin  kd„ 


ln+l  A 

I  Sn+1-Im  sin  k<£n+l"*)  d* 


sin  kd 


n+1 


The  axial  test  dipole  m  is  a  filamentary  electric  line  source. 

In  many  cases  the  surface  dipolo  n  can  also  be  regarded  as  a  filamentary 
V  dipole  located  on  the  composite  cylindrical  surface  S.  In  this 
event  Eq.  19  applies,  and  the  elements  in  our  mutual-impedance  matrix 
are  given  by  classical  induced  emf  theory.  The  field  E*  is  obtained 
from  Eqs.  7-9,  and  the  integrations  in  Eq.  19  are  usually  performed 
by  Gaussian  quadrature  on  a  digital  computer.  When  dipoles  m  and  n 
are  close  together,  however,  a  closed-form  integration  (in  terms  of 
exponential  integrals!  is  employed  for  Z^. 

In  using  Eq.  19  it  is  noted  that  the  numerical  results  depend 
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Co  some  extent  on  the  circumferential  position  of  dipole  n  on 
surface  S.  This  dependence  is  greatest  when  the  angle  of  intersection 
of  segments  n  and  n+1  is  acute.  When  the  intersection  angle  is 
small  it  may  be  necessary  to  average  the  values  obtained  for  Z^  for 
several  circumferential  positions  of  the  surface  dipole  n.  Strictly 
speaking,  should  be  regarded  as  the  mutual  impedance  between  a 
filamentary  axial  dipole  m  and  a  tubular  surface-current  dipole  n. 

Letting  N  represent  the  number  of  unknown  current  samples,  we 
enforce  N  reaction  tests  by  letting  m  =  1,2, •••N  in  Eq.  18.  This 
yields  a  system  of  simultaneous  linear  equations.  This  system  is 
solved  numerically  to  determine  the  current  samples  In.  The 
impedance  matrix  2^nn  is,  of  course,  symmetric.  In  most  problems 
the  antenna  is  driven  by  a  single  one-volt  generator  (Vg*=l)  located  at 
£g,  and  the  antenna  impedance  is  Z  =  1/Ig. 

After  the  current  samples  In  have  been  calculated,  the  far-zone 
antenna  field  patterns  can  be  obtained  in  a  straightforward  manner 
as  shown  in  the  next  section. 


IV.  THE  FAR-ZONE  FIELD  ANALYSIS 

Wt  shall  develop  expressions  for  the  far-zone  field  of  one 
straight  segment  of  the  antenna.  By  superposition,  the  antenna  field 
can  then  be  obtained  by  summing  the  contributions  from  all  the  segments. 

Let  (xj ,y, ,Zj)and  (x2,y2,z2)  denote  the  coordinates  of  the 
endpoints  of  the  segment.  The  current  distribution  on  the  segment  is 
represented  by  1(a)  =  £ I(i)  where  £  is  a  unit  vector  from  (xj.yj.Zj) 

toward  (x2,y2,z2)  and  £  is  the  distance  along  the  segment  measured 
from  the  point  (xi.yj.zj).  In  the  far  zone,  the  vector  potential 
is 


(20)  A  (r,e,<j>) 


tJe~jkr  f6 IU)  eJkfd£ 
4rr  o 


where 


(21)  f  *  x  sin6  cos  $  +  y  sin  0  sin  $  +  z  cos0 


and  d  is  the  length  of  the  segment.  The  source  coordinates 

are  (x,y,z).  The  spherical  coordinates  (r ,  0 , <p)  for  the  observation 

point  are  defined  by  the  following  transformation: 
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(22)  x  ■  r  sin  0  cos  $ 

(23)  y  *  r  sin  0  sin  4> 

(24)  z  ■  r  cos  0 

The  electric  field  intensity  is  E  ■  -  ju>A,  where  the  radial  component 
E  is  to  be  suppressed.  The  rectangular  coordinates  of  any  point  on 
tne  source  are 

(25)  x  *  x\  +  4  cos  o 

(26)  y  ■  yj  +  4  cos  0 

(27)  z  ■  Zj  +  4  cos  y 

where  (cos  a,  cos  6,  cosy)  are  the  direction  cosines  of  the  4  axis. 
From  Eqs.  20  and  25-27 , 

(28)  E«(r,e,<f>)  -  -  jam  e“jkr  ejkfij  i(4)  e;i8Ad4 

4itr  0 

where 

(29)  fj  «*  X)  ;; in  (’  cos  <J>  4-  yj  sin  0  sin  d>  +  zj  cos  0 

(30)  g  *  k  ("cos  .t  sin  0  cos  <t>  +  cos  0  sin  0  sin  <p  +  cos  y  cos  0). 
Noting  that 

A  *  A  * 

(31)  4  *»  x  cos  a  +  y  cos  8  +  z  cos  y 

A  A  A  A 

(32)  0  ■  x  cos  0  cos  $  +  y  cos  0  sin  <f>  -z  sin  0 
and 

(33)  $  *  -x  sin  <}•  +  y  cos  <j>, 
we  find  that 

(34)  E0  *  (cos  a  cos  0  cos  $  +•  cos  0  cos  0  sin  <{>  -cos  y  sin  0)EJi 
and 

(35)  E.  ■  (-cos  a  sin  $  +  cos  8  cos  <J>)E0 

$  X 

Using  the  piecewise-sinusoidal  expansion,  the  current  distri- 
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bution  on  the  segment  Is  given  by 


(36)  I(i)  -  _i -  [I  sin  k(d-i)  +  I,  sin  k£] 

sin  kd  1 

where  Ij  -  1(0)  and  I2  ■  1(d).  From  Eqs.  28  and  36, 


(37)  E  (r,0,*)  -  - 30  .f jkr -  (C  eJkfl 

*■  (1  -  g2/k2)  r  sin  kd 


C*  I2  ejkf2) 


where 


(38)  C  »  sin  gd  -  g  sin  kd  +  j(cos  kd  -  cos  gd) 

(39)  f2  ■  x2  sin  0  cos  $  +  y2  sin  0  sin  +  z2  cos  0 

After  calculating  from  Eq.  37,  the  theta-  and  phi-polarized 
components  of  the  far-zone  field  are  obtained  from  Eqs.  34  and  35. 

After  summing  the  field  contributions  from  all  the  segments 
of  the  antenna,  one  has  the  far-zone  fields  EQ  (r,0,<{>)  and  E  (r,0,<J>). 
Finally,  the  power  gain  is  determined  as  follows  *** 


(40)  Gp(0,*)  - 


(  |e9I2  +  1EJ2  >  r2 

30  |lg|2  R 


where  I  is  the  terminal  current  and  R  is  the  real  part  of  the 
antenna°impedance  Z.  The  directive  gain  is  obtained  from  Eq.  40 
by  replacing  R  with  R,..  The  radiation  resistance  Rr  is  discussed 
In  section  VII.  As  shown  in  the  next  section,  the  effects  of 
lumped  impedance  loading  are  readily  analyzed  with  the  reaction 
formulation. 


V.  LUMPED  LOADS 

Suppose  a  lumped  load  impedance  Zm  is  inserted  in  the  wire 
antenna  at  each  of  the  current  sampling  points  &m.  These  impedances 
may  be  active  or  passive,  identical  or  assorted.  The  current  through 
the  load  Zm  is  denoted  by  Im,  and  the  total  voltage  acting  at  the 
point  £m  is 
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(41)  V.  -  V* 

®  B 


where  V*  represents  a  voltaga  genarator  that  may  ba  locatad  at  tB 
in  aarlfa  with  tha  load  Zq.  From  Eqa.  18  and  41,  tha  ays tarn  of 
linaar  aquations  is 

!,  !n  zm  '  Vi  '  ^ 

n«l 


Transposing  the  last  term,  we  find 

<«)  !  I„  Z»  -V* 

n*l 

Except  for  the  diagonal  elements,  the  new  Impedance  matrix  is  the 
same  as  the  original:  .  The  diagonal  elements  are  given  by 


(44) 


V 

Z„  , 

mm 


mm 


Thus,  the  effect  of  lumped  loading  is  accounted  for  simp  r  by 
adding  the  load  impedances  Zm  to  the  corresponding  diago  A  elements 
in  the  impedance  matrix. 


The  effects  of  the  finite  conductivity  of  the  ante:  na  wire  are 
considered  in  the  next  section. 


VI.  FINITE  CONDUCTIVITY 

After  Ec.  11,  we  assumed  for  simplicity  that  the  wire  antenna 
has  perfect  or  infinite  conductivity.  However,  it  is  not  difficult  to 
account  for  the  finite  conductivity  o  of  the  metal  actually  used  in 
constructing  a  practical  antenna.  Although  the  tangential  component  of 
E  vanishes  at  the  surface  of  a  perfectly  conducting  wire,  it  does  not 
vanish  with  finite  conductivity.  Instead, 


(45)  Et  -  Z8  H* 


where  the  electric  and  magnetic  field  Intensities  and  H* 
are  evaluated  at  the  surface  of  the  wire  and  Z.  denotes  the  surface 
impedance  of  the  wire.  E^  represents  the  longitudinal  electric 
field  component  (parallel  with  the  axis  of  the  wire) ,  and  H±  is  the 

circumferential  magnetic  field  component  (encircling  the  wire) . 

*  » 
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If  each  segment  of  Che  antenna  has  the  same  radius  and  conductivity, 
ve  nay  consider  the  surface  impedance  to  be  a  constant. 

Pro*  Eq.  14  and  Jg '■  n  *  H, 

(46)  H4  -  Js  -  liil 

T  2*a 

where  I(£)  is  the  total  (conduction  and  displacement)  current.  From 
Eq.  45  and  Mg  «*  JJ  *  n, 

(47)  It  -  fa,  -  • 

2"a 

From  Eq.  llf 

N  , 

C«8)  -// 8  Jg-l’ds  +  Us  4‘H’  ds  -  l  Vn  I  (V 

n-1 

From  Eqs.  14,  47  and  48, 

N 

(49)  -ft  1(0  -E  ’  d£  +  Zs  ft  I(Oh’  d£  -  £  V„  l’(£n) 

0  0  w  n»l 

I  f 

In  Eq.  49,  it  is  to  be  noted  that  E^  and  H  represent  the  free-space 
fields  of  the  axial  test  source,  evaluated  on  the  surface  S.  If 
the  wire  radius  is  small,  Eq.  13  yields 

(50)  H*  -  U&L 

♦  2na 

From  Eqs.  49  and  50, 

(51)  -/J  I(0\E'd£  +  Z8  ft  I(£)l’(£)d£  -  l  Vn  It(£n) 

2wa  n"1 

In  reaction  test  m  the  current  function  I*(£)  is  given  by  Eq.  16, 
and  Eq.  51  reduces  to 

(52)  -ft  KO-E*  +  K  I1™’1  I(£)l’  (Od£  -  Vffl 

0“  Vi 
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Using  Eqs.  10  and  16,  we  readily  evaluate  the  last  Integral  in 
Eq.  52  to  obtain  the  sinultaneous  linear  equations: 

N 

<»>  I  Vm  ■  \ 

n»l 


where 


(54)  Z, 


mm 


.  Z__  +  T2  kdm 
mm  8..,.n  - El 

sii 


-  sin  2  kdm  +  2  Vd^-sin  2  kd^ 

rr~. - -  - r-2: 


sin  kd, 


m 


sin  kd^i 


<55>  Zm,®+1  -  Zm>(n+)  4-Zs^  sin  kd^  -kdm4,,  r.,s  kd^ 


sin  kd, 


m+1 


(56)  Zm.m-1  *  zm  r-1  +  & _  sin  kdm  -  kdm  ‘*os  kd 

4*ka  2. 

sin  kd 


m 


'm 


The  mutual  impedances  for  a  perfectly  conducting  antenna  are 
given  by  Eq.  19.  "qo.  54-56  show  the  adjustments  that  must  be  made 
to  obtain  the  mutual  impedance  matrix  2mn  for  the  antenna  with  finite 
conductivity.  Only  the  diagonal  elements  ?mn  and  their  next-door 
neighbors  and  Zm  p,_[  are  affected,  l'or  the  remaining 

elements,  Z  ’  •»  /.  . 

mu  r  n 

It  has  been  pointed  out  that  the  matrix  Zmn  is  symmetric. 

From  Eqs.  54-56  it  can  easily  be  shown  that  is  also  symmetric: 

^mn  "  Znm* 

Previous  investigations  of  finitely  conducting  wires  have 
used  comparatively  crude  methods.  In  some  cases  the  distributed 
effects  of  the  conductivity  have  been  approximated  with  equivalent 
lumped  loads.  In  a  widely  used  perturbation  method,  the  current 
distribution  is  assumed  to  be  the  same  as  that  on  a  perfectly  conducting 
antenna.  Since  our  new  analysis  avoids  the  usual  approximations,  it 
is  believed  to  be  accurate  for  a  wider  range  of  conductivities. 


VII.  RADIATION  EFFICIENCY 

If  the  antenna  has  only  one  voltage  generator  Vg,  the  time- 
average  power  input  is 

(57)  Pin  -  Ug  |2R 
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where  R  is  ths  rssl  pert  of  ehs  ante  «n*  imped  tries  *e  the  terminal 
point  ig.  Ths  rsdistsd  powsr  is 

(SB)  Pr  -  |lg|2Rr 

where  Rr  is  ths  radiation  resistance.  The  power  loss  is 

(59)  rt  -  |ig|\ 

where  Is  the  lo:?s  raalstanca.  From  conservation  o i  enar  y, 

(60)  R  -  \  +  R* 


The  radiation  eiiiciency  Is  defined  by 


(61)  Pr/pln 


Rr  -  R-R* 

tt  nr- 


With  finite,  conductivity  and  lumped  loads,  the  antenns  impedance 
(and  the  resistance  R)  can  be  calculated  with  the  techniques  developed 
in  the  preceding  sections.  It  is  assumed  that  the  antenna  is 
situated  in  free  space  or  located  over  a  perfectly  conducting  ground 
plane.  Therefore,  power  loss  occurs  only  In  the  lumped  loads  and  the 
antenna  wire  Itself.  The  power  loss  in  the  lumped  loads  is  given  by 


(62)  PM 


N 

■  l 

m-1 


llj  Rm 


where  is  the  real  part  of  the  load  impedance  Zm. 

Now  let  us  consider  the  power  Pd  dissipated  in  the  antenna  wire 
with  finite  conductivity.  If  E  and  H  denote  the  rms  field  inten¬ 
sities,  the  time-average  dissipated  power  is 


(63)  Pd  -  -  Real  #  <E  x  H*)-ds 

0  mmm 


where  S  denotes  the  antenna  surface.  In  view  of  Eqs.  45  and  46,  Eq.  63 
reduces  to 
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(&*)  *<j  *  ?2_  fL  1  1*  dt 

2ira  0 

where  the  surface  resistance  Rg  is  the  real  part  of  Z8.  Eq.  64  can 
be  written  as  follows: 

N+l  2n 

(65)  Pd  -  Rs  J  f  1(1)  1*  U)  d£ 

lira  n-1 

where  the  current  on  segment  n  is 

(66)  I(£)  ■  ln_i  sin  k(£n-fc)  +  In  sin  k(&-£n-1) 

sin  k  dn 

From  Eqs.  65  and  66,  the  dissipated  power  is 
N+l  r 

(67)  rd-?s  ?  — ir-  (I„.1i;.1+lnO(2kdn-8ln2kV 

8irka  n=l  sin  kdnL 

+  «Vl  +  Cl^Xsi"  kdn-kdnC08  kd.7| 

The  total  power  lois  '  s 


(68)  P£  -  P££  +  Pd 

From  Eq.  59,  the  loss  resistance  is 

(69)  R,  -  Pt/  |lg|2 

Finally,  the  radiation  efficiency  is  calculated  as  shown  in  Eq.  61. 

An  alternative  method  is  to  determine  the  radiated  power  by 
integrating  the  power  density  in  the  far-zone  field.  However,  the 
technique  developed  above  is  much  faster. 

When  R  and  R^have  been  computed,  the  radiation  resistance  Rr 
is  obtained  from  Eq.  60.  The  directive  gain  Gd  of  the  antenna  can 
then  be  calculated  from  Eq.  40  with  R  replaced  by  Rr. 
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VIZI.  COMPARISON  OF  CALCULATED  AND  MEASURED  RESULTS 


Th«  Smith  charts  in  Figures  3,  4  and  5  compare  tha  calculated  and 
measured  impedances  for  three  different  antenna  configurations.  Although 
these  antennas  are  reasonably  simple,  no  adequate  theory  or  com¬ 
putational  method  has  heretofore  been  available. 

The  entennas  were  constructed  of  brass  rod  with  a  diameter  o 
1/8  inch.  Each  antenna  has  straight  segments  with  a  right-angle 
bend  every  5  cm.  The  antenna  was  mounted  at  the  center  of  a  large 
aluminum  ground  plane  and  fed  with  a  50-ohm  coaxial  cable.  A 
General  Radio  UHF  Admittance  Meter  was  used  for  the  measurements. 


It  ms/  be  noted  in  Figures  3-5  that  the  calculated  and  measured 
impedances  show  good  agreement. 

Figure  6  shows  the  Impedance  versus  frequency  for  a  16-sided 
polygon-loop  antenna.  These  data  were  calculated  with  the  technique 
described  earlier  in  this  report,  except  that  the  current  function 
1(1  )  is  not  assumed  to  vanish  at  any  points  on  a  loop  antenna.  F  >r 
comparison.  Figure  6  also  shows  Storer's^  calculations  for  a  clrcu  ar 
loop  with  a  circumference  equal  to  the  perimeter  of  the  polygon- 
loop.  The  close  agreement  between  our  calculations  and  Storer's 
is  a  further  indication  of  the  accuracy  of  the  new  formulation. 


Figure  7  shows  measured  and  calculated  impedance  versus 
frequency  for  a  square  halo  antenna  driven  from  a  50  ohm  coaxial 
cable  through  a  ground  plane.  The  good  agreement  indicates  that 
our  technique  is  accurate  for  three-dimensional  loops  as  well  as 
dipoles . 

Figure  8  shows  the  radiation  efficiency  versus  frequency  for 
the  antenna  in  Figure  4,  calculated  with  the  techniques  developed  in 
Section  VI.  The  antenna  in  Figure  3  has  a  slightly  higher 
efficiency,  and  the  antenna  in  Figure  5  has  a  slightly  lower 
efficiency. 


IX.  CONCLUSIONS 

The  experimental  approach  to  antenna  measurement  and  optimization 
is  a  time-consuming  cut-and-try  affair  even  with  modern  swept- 
frequtney  equipment.  To  augment  the  experimental  approach,  we  have 
developed  an  accurate  and  efficient  digital  computer  program  for  e 
wide  variety  of  wire  antenna  configurations.  This  program  provides 
a  complete  analysis  of  the  antenna  characteristics  including  the 
Impedance,  current  distribution,  radiation  efficiency,  field 
patterns  and  gain. 
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Fig.  7.  Impedance  vs.  frequency  for  square  halo  antenna. 
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RADIATION  EFFICIENCY  (PERCENT) 


Fig.  8.  Radiation  efficiency  vs.  frequency  for  antenna 
in  Fig.  4. 


This  report  develops  the  theoretical  basis  for  our  new  for¬ 
mulation,  including  the  effects  of  multiple  excitation,  lumped 
loading  and  finite  conductivity.  Experimental  impedance  data  are 
included  for  several  three-dimensional  antennas  to  verify  the  accuracy 
of  the  calculations.  Appendices  present  the  computer  program  and 
an  analysis  of  scattering  by  three-dimensional  thin-wire  antennas. 

The  computer  program  is  convenient  to  use.  One  merely  gives 
input  data  to  specify  the  antenna  shape,  wire  radius,  con¬ 
ductivity,  frequency,  lumped  loads  and  generators. 

The  program  is  quite  efficient.  With  an  IBM  7094  computer,  an 
analysis  of  the  antenna  in  Figure  3  requires  only  0.5  seconds  for 
each  frequency.  The  antennas  in  Figures  4  and  5  require  1.0  and  1.5 
seconds  respectively. 

We  have  found  that  these  techniques  for  wire  antennas  in  free 
space  can  be  applied  also  to  a  more  advanced  problem  of  considerable 
curren.  interest:  A  wire  antenna  exciting  a  circumferential  slot 
on  a  rocket  or  projectile.  Continuing  investigation  in  this  direction 
is  recommended. 


ACKNOWLEDGEMENT" 

The  experimental  results  reported  herein  on  antenna  impedance 
were  measured  by  Danai  Lekhyananda  and  George  A.  Richards.  Their 
assistance  is  sj.i.  <  roly  appreciated. 


125 


REFERENCES 


1.  Schelkunoff,  S.  A.  and  Friis,  H.  T.,  Antennas,  Theory  and 
Practice .  New  York:  Wiley,  1952,  p.  401. 

2.  Schelkunoff,  S.  A.,  "On  Diffraction  and  Radiation  of  Electro¬ 
magnetic  W a ves",  Physical  Review,  Vol.  56,  August  15,  1939. 

3.  Carson,  J.  R. ,  "Reciprocal  Theorems  in  Radio  Communication", 
IRE.  Proc.,  Vol,  17,  June  1929,  p.  952. 

4.  Storer,  J.  E.,  "Impedance  of  Thin-Wire  Loop  Antennas", 

Trans.  A.I.E.E.,  Part  I,  Vol.  75,  Nov.  1956,  pp.  606-619. 


126 


****»?*.+ 


APPENDIX  I 
Scattering 

In  the  scattering  problem,  we  assume  the  antenna  has.  perfect 
conductivity,  no  voltage  generators  and  no  lumped  loads.  Let  the 
Incident  field  (E^,Hj.)  be  generated  by  an  electric  current  source 
J-q .  Then  ■  0  and  Eq.  11  reduces  to 

<*»  -!!s  4 -s’  <i’  -///„  4  '  i'  iv 

As  in  Section  111,  Eq.  70  leads  to  a  system  of  simultaneous 
linear  equations: 

N 

(71)  l  lnZm‘  ///„  Ji-E,  dv  -  vm 
n*l 

where  the  mutual  impedances  Zmn  are  given  by  Eq.  19. 

The  simplest  case,  and  one  of  considerable  interest,  is  that 
in  which  the  primary  source  is  at  a  great  distance  r  from  the  antenna. 
In  this  case  the  field  E*  of  test  dipole  m  is  given  by  Eqs.  34,35 
and  37  in  Section  IV.  ~IL 

To  analyze  the  plane-wave  scattering  properties  of  a  three- 
dimensional  thin-wire  antenna,  we  let  the  source  Jj_  be  remote  and 
infinitesinal  in  extent.  Then  the  right-hand  side  of  Eq.  71  reduces 
to 


02)  V,  -  i  •  JJJV  Jjdv  -  S’  •  4 

where  the  incident  field  Ej_  is  evaluated  at  the  coordinate  origin 
which  is  at  or  near  the  wire  antenna.  To  illuminate  the  antenna 
with  a  plane  wave  arriving  from  direction  ,fwe  locate  the 

source  JL  at  the  position  (r^,^)  and  evaluate  ^  in  Eq.  72  at 
(r,0j.  .  Likewise,  we  orient  E^  in  the  9  or  £  direction  to 

investigate  the  scattering  for  a  theta-polarized  or  a  phi-polarized 
incident  plane  wave. 

It  may  be  noted  from  Eqs.  71  and  72  that  the  solution  for  the 
scattering  problem  proceeds  in  precisely  the  same  manner  as  in  the 
antenna  radiation  problem.  Indeed,  the  currents  In  induced  on  the 
antenna  in  the  scattering  problem  are  the  same  as  those  excited  by 
a  series  of  equivalent  voltage  generators  Vm  inserted  at  the  current 
sampling  points  £m.  The  voltages  of  these  equivalent  generators 
are  given  by  Eq.  72. 
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When  the  solution  has  been  completed  for  the  simultaneous 
equations  (Eq.  71),  Eqs.  34,  35  and  37  may  be  employed  to  calculate 
the  scattered  field  Eg  at  any  distant  observation  point  (r , 8S , 4>s ) . 
The  monostatic  or  bistatic  echo  area  of  the  wire  is  then  given  by 


(73)  o  -  Aim 

r*» 


If  the  wire  has  lumped  loads  and  finite  conductivity,  the 
effects  are  readily  treated  in  the  scattering  problem  with  the 
techniques  developed  in  Sections  V  and  VI. 


L 


E<-E* 
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APPENDIX  II 
The  Computer  Program 

Figure  9  presents  the  main  computer  program  (in  Fortran  IV) 
for  three-dimensional  dipole  antennas  in  free  space.  A  similar 
program  is  available  for  three-dimensional  loop  antennas. 

The  symbols  for  the  input  data  are  defined  as  follows: 

IND  1  if  patterns  are  desired,  or  0  if  not  desired. 

NCAS  Number  of  frequencies  to  be  analyzed. 

NG  Index  number  of  antenna  driving-point  terminal. 

NP  Number  of  points  on  the  antenna. 

ACM  Outer  radius  of  wire  in  centimeters 

BCM  Inner  radius  of  wire  in  centimeters 

CMM  Conductivity  of  wire  in  megamhos/m 

(-1.  for  perfect  conductivity). 

EPSIR  Relative  permittivity  of  the  internal  medium. 

FMC  Initial  frequency  is  MHz. 

DFMC  Frequency  increment  in  MHz. 

XC(I), 

YC(I), 

ZC(I)  The  coordinates  in  centimeters  of  point  i  on  the 
antenna  axis. 

VG(I)  The  complex  voltage  of  the  generator  Vj\ 

ZLD(I)  The  complex  lumped  load  impedance 

Typical  input  data  are  shown  in  Fig.  8  at  the  end  of  the 
program  listing.  These  data  apply  to  the  antenna  shown  in  Fig.  3. 

In  this  example,  far-zone  pattern  calculations  are  not  desired. 

The  impedance  is  to  be  calculated  at  19  frequencies  starting  at 
IOC  MHz  with  increments  of  50  MHz.  The  antenna  has  no  lumped 
loads.  There  is  just  one  voltage  generator  and  it  is  located  at 
the  center  of  the  antenna.  The  antenna  shape  is  defined  by  an 
ordered  list  of  the  coordinates  (x^,y^,z^)  of  five  points 
(NP*5)  on  the  antenna  axis,  where  i  *  1,2, 3,4,5.  There  are  three 
current  samples  I l ,  I2  and  I3  to  be  calculated.  The  unit  voltage 
generator  is  located  at  the  point  where  the  current  I2  is  defined,  and 
therefore,  NG  ■  2.  The  wire  radius  is  0.1588  cm  and  the  conductivity 
is  10.4  meganhos  per  meter.  Since  the  wire  is  solid,  the  inner 
radius  BCM  is  zero. 

In  the  program,  WAVE  represents  the  wavelength  in  centimeters, 

AK  *  ka  and  TPL  **  k.  The  surface  impedance  ZS  is  calculated  in  the 
subroutine  ZHOLOW^  Subroutine  TDD  calculates  the  mutual  impedances 
Zjpn  (Eq.  19)  and  Zmn  (Eqs.  54-56).  In  the  program  Z  is  denoted 
by  C(I,J).  Just  above  statement  100  the  lumped  loads  are  accounted 
for  in  accordance  with  Eq.  44,  and  the  impedance  matrix  is  made 
symmetric. 
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SEXELUTE  IHJOR 

SIHJOO  GO 

SIBFTc  MAIN  NOOECK 

1  FORMA T  (  1  X *~it-  1  5*7} 

2  FORMAT  <  J  X  » tip  I  B»  7 1 

-  5  FORMAT 1 1  HO  I 

|  7  FORMAT (7F10.5) 

\  8  FORMAT  <7 1 1  i<  I 

|  It  FORMAT  UH|) 

5  I?  FORMAT <1X*HF|S*1) 

COMPLEX  CN«CL 

i  COMPLEX  CJ..7Lr>(50)«VG<50I 

I  COMPLEX  ET1 *rT?.SPt*EP2 

i  COMPLEX  h  T“>«*  P*-I*C  TT  <5  'I  «EPP(  50  ) 

COMPLEX  Z1  !*ZI?*VI!*Y12 
v  COMPLEX  C. 

i  COMPLEX  0'11*CP!2«CP21«CP22 

>  COMPLEX  2HCLOW*7S«ZH 

DIMENSION  X<Bn>.Y<'>")«ZI5>0 

DIMENSION  :>IC!0)«VC50)«WC5'i)«CPS<,SG>«SPSI50> 
DIMENSION  XC (50)* YClb  1) *ZC I5G ) 

DIMENSION  CUK  <50  )  «  SDK  <59  )*5T!>K<50| 

TP**6«2H3 I  M->3 

READI5.0 > |ND*NC«S«NG.NP 

WRITE (6*8 »iND«NG.NP 

WRITE<6«5) 

READ IS* 7  |ACM*eCM»CMM«FPStR«FMC*DFMC 
WRITE (6*2 )ACMtBCM«CMM«ePSIRtKMC 
WRI  TE  (6*  5  ) 

NM*NP- I 

N»NP-2 

DO  20  I « I *NP 

REA0IS.7IXC ( | )«YC( I )*2C( I  I 

F I  •  I 

20  WRITE(6*2|F|«XC(I)*YC(I  )«ZCII) 

WR I T£ (6«S ) 

DO  25  |«I#N 
READ<5<7IVG<I )*2L0< I » 

F I  *  I 

25  WPITEI6«?)r I • VG< I) *2LD( I > 

WRITE (6«5 I 

00  *00  NCASt«l ,NCAS 

WAVE»30000«/FMC 

AL*ACM/WAVC 

ak«tp*al 

tpl*tp/wave 

DO  28  l«!«NP 
X ( I l«XC< I l*TPL 
Y <  I  »«YC<  I  >*TP|_ 

28  Z < I |»ZC  < I  )*TPl 
aL*ncM/wAv: 

HK*TP*8L 
FCPS«FM'-*|.E6 
SIGMAsCMM#!  *56 

ZS»l*0«.H) 

IFICMM.c.T*  ‘  •  )Z'>«7h0L0v(A<*mK«FCPS»SI6MA*FPS!W) 
WRI TE  < A*P )/* 

W«!Tr<6*5) 

ZH«ZS/«*/TP/Al< 

call  TO.)<A<*C*Cf"-1,n*X«  Y«/*ZH«NV ) 

DO  100  I*| «N 
Cll«l  )*F<  |  *  |  ) 4 7L>' I  I  I 
00  IO0  J*!*N 

F|«| 


r*  n  n 


p9 


rj«j 

CIJ.I  )■<*(«  «J> 

100  CONTINUE 
WR I TE (6*5  ) 

CROUT  METHOD  FOR  SYMMETRIC  matrix 

N  REPRESENTS  THE  NUMtiER  OF  SIMULTANEOUS  LINEAR  EQUATIONS 
SET  UP  THE  SQUARE  MATRIX  C(I<J)  BEFORE  ENTERING  CROUT  PROGRAM 
NN«N*I 

DO  120  L«1«N 
ILL-L-l 
00  120  I *L  »N 

if<l«eo*i igo  to  lie 

00  I  17  K* 1 *LLL 

117  C  (L  .  I  >  *0 (L  *  I  )-C(L«K)*C(K.  1  I 
lie  CONTINUE 

IFlL»t0.1  K-0  TO  120 

c ( I  «L)*r  ii.«  n 
C  (L  «  I  >««“<!_•  I  I/CIL.LJ 
120  CONTlNuf 

C  AT  THIS  POINT  CALCULATE  OR  READ  THE  RIGHT-HAND  COLUMN  C(I.NN) 

00  124  1  =  1  *N 
124  C  < I « NN ) *  V  0 (  I  t 
DO  130  L=1 «N 
l LL*L- I 

IF(L.E0.  I  ).*,0  TO  120 
00  127  k  ••  1  « L L L 

1 27  C  <L  «NN )  =C  (L  «NN  )-C(L#K)*C  (X  *NN  ) 

128  CONTINUE 

130  C(L«NN)=..  <L.MN)/C(l.«L> 

IFIN.E'O.  1  )G0  TO  138 
00  136  L«2»N 
1 *NN-L 
I  1 » 1  +  1 

00  136  K  =  I  I ♦ N 

136  Cl  I .NN)=C ( I «NN)-C(  I *K>*CtKlNN) 

138  CONTINUE 

C  NOW  THF  SOLUTION  IS  STORED  IN  THE  RIGHT-HAND  COLUMN  C(I*NN) 

00  140  1  =  1. N 
CJ(l)=C(I.NN) 

CMAG«CAHS(CJI I  )  ) 

CRCAL  =  f  EALICJI I  ) > 

C1MAG=AIMAG<CJ< I > 1 
PHASE=.0 

IFICMAG.EO.O. >G0  TO  139 
PHASE  =  S7 .29576* A T AN2 I C I  MAG .CRFAL  > 

139  WR I TE  <  6. 2 ) CM AG .PHASE 

140  CONTINUE 
WR I TE  < 6.5  ) 

Z1 1*1  ./TJ(NG) 

R*RCAl(Z1 1 > 

RL».0 

IF (CMM.lT .0. )G0  TO  1 60 
PD-.O 

00  150  1=1. NM 
0K=0< I ) 

SOK*SlN(OK  1 
CDK*C0S(0K ) 

ST0K«SIN(?.-*DK> 

CN«  (.0.  .v) ) 

CL*  (  .0.  ••'> ) 

I F ( I.GT.I >CL*C J ( I - 1 > 

I F ( l.LT.NM)CN*CJ(  I  1 

PA=  (2.*OF-!.TD<)*(CA-*GICL)*«?+CAHS(CNl**2) 

PH  =  4.*  <  SO<-OK*CnK  )•  <PEAL<CL)*Pf- AL(CNl  +  AlMA5(fL)*AIMAG(CN)  ) 
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150  P0»P0'MPA+PB}/$0K/r»0K 
RS*REAL<ZS> 

PD*RS*PD/ I 4 • *TR*A< ) 

RL*PO/trArtS(C«MNG) I  »*?> 

160  EFF 1C* 1 00*  *  <R-PL ) /U 

WR l TE (6*  2  )FMC«KFF ! C  tZl 1 
IF ( 1NO*EO«0 )G0  TO  350 
WRITE (6* 1 1  ) 

CJG*CJ(NG) 

00  206  I  ■  1  «N 
206  CJ( t 1*CJ< I l/CJG 
00  340  NPP»1«3 
TH«90« 

PM*#0 

1F<NPR«EO.T)oh«RO. 

00  300  NT  * I  *  19 
THR»«01745'29*TH 
PMR«#01745.52o«ph 
CTM«COS<THm  j 
STH«SIN(THR> 

CPM«COS<PnH) 

SPM«S|N(PHR) 

DO  221  1*1 «N 
ETT ( I >*< .0*«0  > 

EPP( I >*(.0«.0) 

221  CONTINUE 

00  22 5  K*1 «NM 
XA«XOC) 

YA*Y<K» 

ZA*Z<K) 

KP*K«-1 

Xa*X(KP) 

VB*Y<KO> 

ZB-ZOCPJ 
DD*D(K 1 

CALL  ZFF (XA« YA« /A* Ybt ZB  *  DO* CTH* STH« CPM« SPH* ET1 «ET2»EP| «FP2> 
ETT<K)*ETTIK)4FT2 

epp  <*  i  «epp  <*■ )  4rDr> 

IFOC.EO.l >G0  TO  ?2b 
KM*K>1 

ETT<KM1«ETT(K'11*ET1 
FPP(KM)trpp(<M>4.EPi 
225  CONTINUE 

ETH*  <  »0«  »0 1 
£PM* ( »0« «0 ) 

DO  250  1*1 *N 
ETH*ETH*CJ<  I  )*FTTU  ) 

EPM«EPH4CJ<! >#FPP<I » 

250  CONTINUE 

OPT*  (CAbr><r.TH>«*2)/R/30« 

CPP*  (CABS  <£PH  J  *t*2  >/R/30» 

GP*GPT4GPP 

WRI TE (6«  2  1TH<PM<GP<GPT <GPP 
IF  <NPR»EO< 1 )PH«PM*20. 

I F  < NPR <  E  Q  <  1 >00  TO  30o 
TH«TH-|0. 

I F  «NT»lT •  1  >'  )G0  ?0  300 
TH* 1 0* (NT -9 ) 

PM=1B0< 

IF<NPP<KQ< 3 IPMs-90, 

300  CONTINUE 
340  WRITE<<S<M 
3bO  CONTINUE 

WRITE<fc<M 
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IPJ  IND.NE.  •  >' )  w>J |  T£  (6. 1 1  ) 

T I ME “CLOCK  <0#  0 1 6667  > 

WB!TE<6.?)T|MC 

400  PMC«FMC+f>EM> 

STOP 

END 
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The  simultaneous  linear  equations  (Eq.  18)  are  solved  between 
statements  100  and  138.  The  voltage  generators  are  accounted  for 
in  statement  124.  Statement  3  3**  writes  out  the  magnitudes  and 
phases  of  the  complex  current  samples  I.  »  CJ(I).  The  dissipated 
power  is  calculated  between  statemencs  140  and  160  using  Eq.  67. 
Iheradiaticn  efficiency  EFFIC  calculated  in  statement  160 
neglects  the  power  loss  in  the  lumped  loads  and  assumes  there 
is  just  one  voltage  generator.  The  antenna  impedance  Zll  is 
calculated  under  the  assumption  that  there  is  just  one  generator 
and  that  it  has  unit  voltage.  The  frequency,  efficiency  and 
impedance  are  written  out  after  statement  160. 

The  far-field  patterns  are  calculated  between  statements  206 
and  340.  Specifically,  the  power  gain  G?(Eq.  40)  is  calculated 
in  the  three  principal  planes  with  10-degree  increments  in  6  and 
20-degree  increments  in  $.  In  addition,  the  theta-polarized  and 
phi-polarized  power  patterns  GPT  and  GPP  are  calculated.  These 
calculations  use  subroutine  ZFF  which  is  based  on  Eqs.34,  35 
and  37. 


